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ABSTRACT 


This report is primarily concerned with the 
modification of the mean temperature profile in the surface 
boundary layer of the atmosphere caused by a change in 
surface temperature. 

A review of the literature relevant to this 
subject is given. An experiment is described in which a 
change in surface temperature was created by irrigation of a 
Small area of grass, and its effect on the airflow measured. 

Predictions of an analytical local-advection model 
were eeicuea with the observed temperature-profile 
modification. In individual experiments the comparison was 
inconclusive because of inaccuracy in the experimental data. 
However, data were averaged over all experiments and there 
was found to be no Significant difference between model and 
experimental values of the shape factor used for comparison. 
The model average temperature-gradient stabilisation was 
found to agree reasonably closely with the experimental 
averages. 

Wind-profile modification was observed to be 
stronger than that expected from a numerical local-advection 
model in both several individual experiments and on average. 

Several improvements in experimental technique are 


Suggested for consideration by future experimenters. 


Ly 


oveys zi. OSES 


is pb 


te Pimgx oe) ‘gi \taneinipaend 


[ ah 7 * a 


neko sat et god senoh? 
: ; ‘Le 
| 7 7 


ay 


ACKNOWLEDGEMENTS 


IpWwishetoechank ir. Kod. Hage, for, his patient rand 
helpful supervision of my study, and his help in execution 
of the temperature-profile modification experiment. 

iWwanvgraterul to Dr. J._L. Honsaker, who. designed 
the temperature-sensor network and helped carry out the 
experiment. 

Financial support for the field experiment was 
provided by the National Research Council. Most of the 
meteorological field equipment was provided under grants 
from the National Research Council to the Institute of Earth 
and Pianetary Physics, and by the Water Resources Research 
Program of the Inland Waters Directorate of Environment 
Canada. 

I wish to thank the Department of Soil Science for 
permission to use the Ellerslie site, and Dr. R. B. Charlton 
and Dr. J. A. Toogood for serving on my examining committee. 

Finally I wish to express my appreciation for the 
financial support I have recieved through a Canadian 


Commonwealth Scholarship. 


- 


‘ ¥ 
beavyoxe. og jecw 
7 


[ ; a : eu : _ ian 1 
Woive7) eed 3 \-Saoqqive tak 
, 5 7 ; in ; ae oes we 
i 7 
4 sbdersiod lll o, alata) t 


W ve ti 


TABLE OF CONTENTS 


Page 
ABSTRACT se eecececcccesccecce sce e sce ee sec ees ccc ccccces iv 
PIOKIOW LEDGE MISNT 50 cre a aie o eielaleleie os <'elels 6a e100) 4010s o1s10/4'6/0's's/s om V 
Dim tr aCO NIE NI Sle lelsistew tele sine a esie os «isle ws iasiclciais «619,90 «0 6s vi 
ee ei Canty ACl 171, Date tetel alate! sel ais aie ajels al «sisal aisle ale sie ele) ele sie iclalele «eV ti): 
LIST OF, PIGURE SIRS ac 6 sve sivieiels pcos cle o s cic ss cis ei siscle os ces 2'e XV 
CHAPTER 
if UENTRODUCT ION ( ccicatoia «0.0 sie © 0's sco 0 4 010 sls os 00 <0 0 1 


ita! THEORETICAL BACKGROUND FOR AN INTERNAL 
BOUNDARY LAYER seaeeuseeesee0aea”aeasesesee@eeoeeeseseeteegee@#*e#ensnve2e8ee8 6 


Pes GOV 2 iait in Gand a G01) Cleats cictclaisialcletale)sie\wicis ela a ost! 6 
rere Elimination of Velocity as a Variable ..... 8 


ees Power Law Wind and Diffusivity 
Profiles @eeoeeaoeeoeeeev#weeeeeseseeeoe#seeeeeee oaeseeseeee 9 


LET ANALYTICAL INTERNAL BOUNDARY LAYER 
MODELS @e@aewzeoe*e@esee*ee2nescseseaeveeeaeeoeenekeeesreeee#ee#eeeeeeeseses#*eeaeesee? 1 


Bett Theshogel (Om sSil CO Neher < S o.atewrt sale "Se bes cu 0 ees 12 


3.2 The Microclimate Modification Model of 
de Vries @®aeaoeoeee@eesPcege@eese72#eee#@eeed@e@ese*e#se@*se@e#es#sesteeeeeseseeskee @ 14 


os} The Model of Rider, Philip, and 
Bradley eeee#ss#see?e#*eaees?seenmeenteeeoeseeseee@exn,newsetgeess¢9essecm68t— et ®t @ @ 16 


3.4 Limitations of the Rider, Philip, and 
Bad ev mOUC Te sc wrsisin soe a alalnleteie(sietsinisiels #in iets e's 23 


Vi 


Vs eo 


» a 5 7 
¥ cine bh Ce ee Hee er eee © RB ORS US ere Yee aeeanies i 
; * ae 7 
7 ‘ evn q P's 
cv ~¢ Oe Sebi @eor ss aweere ed erep ree? sv FA . oD. €0 etic 
; Nae *, ' 


ra a 


Jvchlu 4s biveh New Redanesen en vaMEnebe #3 atgean*4@ 


6 ews a4 ooete ee 8a 8 cae Wave tleuae ens acs es QSR00ER See oo 


i 
7 e ts 
‘ Py : ah | 
Hae odie eevee avrs@ os ae @ « “*@ ACTOR HAHA P 7 
LAWRSE TMT, HA BOT. detiortcQde Tea estoy 1 a 


sateen edieves sends See MO eAd THAd4aGE 
{ y A 
cove udaceedenweagmes, BROI@RNDS Pareseves 
i ose @fteltev a BF YFfoeLToV Yo nolsen resis) , 


giivianitid faa bark ws) 3eNoF 


id ++ aoe ar ee were de te pes e pes eee er els asf lyon - 7 
? ote 
$2%K.S YRAOWOOR EMIT CES GdoTTTAAMe 
a 79h 0h SD OS SO Aen ie Dee ee a ee 2.5900 


t ‘somr aaht <emenmpene es oosdat 29 Letok od? 


to Lebon solaoett bot etoaifpoinin aft 
fe] ‘sheen ee ee wee we oh adie Gh akin See ww ew a: ee lay eb. Whe 


0 hea {qhnens \hebta 30 rebom 
ot «hn ay ea ean 


! Teo mag 
oe ’ phic sig 
seat Ser : , wa. 


CHAPTER Page 


Se2 TReBNOd ClLOE PYAR As sre tales 6 Scat < Deedes secs 25 
36 MhOeo Nn WCC ZAMNOEHOGEs sieracte sic clevclels.s Stele sjale.ce «es 26 
ed: Comparison of Shwetz and Philip 
SOU Ur POs im aare = ciciclevel cial calc velele's ies « 0.0/2 00's +10 s'e 29 
iby: NUMERICAL INTERNAL BOUNDARY LAYER 


MODELS @a@eaeeesesstesdesedeeteease@eeeeoewoeaeetenremeseeseaeaeexn2uoeveeseesd=#@ 33 


4. Liem xan Geuengtn Nodel of | Taylor). <<. «=. 33 
4.2 Weisman's Extension of the Taylor 
WOW Talk, lag boo oH Ch Gini AO OO Dee oe 36 
eS} ThesModelvot Rao, Wyndaasda, and Cote ....'... SY, 
4.4 DhemModelBoheNOV1 KOVawmenals 4s sto so o's clos ses sce Sie, 
V EXPERIMENTAL STUDIES OF 
PEALE RREURESPROFELE MODIFICATION. Fowisie wane 4O 
So The Experiment of Rider, Philip, and 
PP BEOQGlCVise sete eiele s'eic.e ss elcle sities olv «is/e/u eles os i= sis 40 
ete The Observations of Dyer and Crawford ..... 43 
ine INSTRUMENTATION, SLE, AND PROCEDURE 
OFF TH EeteGERSOLILESSXPERIMNENT <i ean oe sn cae cee Ky 
6. The Temperature SeENSOLS ..caccccccccencevcecs Ky 
6.2 The Temperature-Sensor Network seccesecccee 46 
6.3 The Sonic Anemometer-Thermometer ...ecceeee 52 
6.4 The Shear Stress Meter and 
Ua) Oe lente tatal alate le ata'a|sicleis!s elolaisisialsleleieic!efsis 53 
O65 Retrieval of Data from Magnetic Tape ...... a6 
ewesls. Bepe te Wet aan ts C Calcrets teal sls aie isi clots als lalalainiais olsiere 56 
6.7 TRSELUMEN EL GAVYOUC. Sine cuisine wccericewncesicas OF 
6.8 Experimental Procedure eeeceeeereescceecece 59 
fame TEMPERATURE~PROFILE MODEFICATION sclsswcicess 63 


Vid 


awe ee eae : a3 ial 


‘voeerdm ae od : bpags 7 > Hovde oF 
4 Bai), 5 7 Nes site ee be et ta ic 


an a i os < uhh “rowers 
7 penn aes 
THAGRCON vaya 


ne ot wales 
, 


’ ey. sag : 
* a 
A 
‘ i] a) , 
. ’ . 4 
1 
Vi ~ 


I % . % 
; 
of ~< 
' a 
’ ® 
LaF i b 6 
de» | : 
ro ® nm es 
* e«s - ea | 
7 


: ae ral bees 
Prreven Me eyed apwnasugs 


’ 7 i) o -_ 


Lz 
huge«da 


le 


a 


CHAPTER Page 
leet Temperature-Profile Modification Data ..... 63 


ere Calculation of Parameters Required by 
the kiden, fr helip, ;and Bradley Model) cscs 82 


Tes Comparison of Observed 
Temperature-Profile Modification with 
the Rider, Philip, and Bradley Model ...... 85 


7.4 Comparison of Model and Experimental 
Stabilisation @@*esteeeea?ee*e@estestsed9s3se@e@eeeesv09seeseecgcgeee@e29eae83228208282 92 
ViGT MODIWIGATION, OBRTHEGVELOCITY , PROFILE 
AND SMUG NET RADIATIVE FLUX ec cicie cc es ae ccs sie 96 
Seat EXPL INetcateWANGHPEOLTENCS sc cccle ewe cneeeiee 96 
De 2 A Simple View of Wind-Profile 
Modification @eev#wre&ee?e@eeaexzs3oe@eaeeseseeseeeneenoeewgegesese2s2neeeasves8seeee 103 
g.3 Modification of the Net Radiative Flux .... 105 
g2).4 CONCLUSION @eseeseeeaeeee?es2se92st@Pte eoeeaseeseeeesesesttse8tmUHDmCUCMAGTCUCMARCUMUCUCMVCUCOCUCB 108 
O51 SUmMAaGymOTe EL NGLAGS. Ys Saye) cate simietsiese Mela s wie cess |) 10S 


Oe2 Recommendations @ea@ameenaee_es#se*esestee?eanae9seteseetre seeoe2esestz 3898 110 
REFERENCES @®@se@ee*e*ee”@*e*ee?eseensnsseeeeoenteoeeeaeaeeeeaeseaeses2s eee de 


APPENDIX {Nomenclature) *eaea@eaaeeee@eeseeeeeseeee#s2#ss&teee2e ® 16 


viil 


2 e 7 a 
s a AR - 
a! iy) ; :| 


in) Be 


' 
: 7 1% 
‘7 ' iy | a. ra : ” 
* a 7 : - i : . me oy 


Vw beak ‘ah ankaaaih ible ee ye 


yd bealopey srercan tat to, nelfretes 
18 a eeeme is on cen ae ve \. 


le lay 
xe 
tikes nokJenka thor oF 
id “oe een @ # febor ¥oTi: sat Tet +4 
i 


to mp: 
eh ee ; o hist E41 ‘ LP 
tetpenbsaqxe tan tenon Ae aap +e rane | 


ae © pie mk © ss) a ee wr Oe Ole oe OP ORES 


«- $eedon8 VSrboasy Swe  ORTRDTES 
at ose be @@ Cikie ee ha ew ut SVITAICAT i 2h | 


; sy 
av ‘e*eeaeeaneen PF eovFr 4 ¥, abiyeess Lert Lie 43 agg coeenit | * deans 


>" ee Be ¥¥-baia Yo. «see faoles ee - 
Eth 4.0486 seeese Oe one leeudmebseice & 616 nH ont 2 7bar .% 


pu, 


i a » “ ‘ ‘ * . , - 
e€ ees THEE) OME Rha? Pen ens ee Rol tha 5 tSEO8 Raa 
BOT rea ae * ? « sedicall afer @e@enetBtGerne vee v?e’s ROTEOAGRO? 


Ol cower dedtivs vedi weresae SEREHECY Gp Po eeeee 


pf * + aieebe nal oie a MSR TREROR ae s.2% 
i S 7 * 
. oe ene €@S6 @e @€ 8 © 46.6. 8») €.¢ @ 2Tonr sone 
——— . ~~ —~— ‘ar. 


33h! Gin ne ee) erreeaen 44c 
- a \ Ny 7 


List OF TABLES 


Page 
Experimental Conditions @eenseeeeaseesee*ees?e#2eeseteeeesee#e 64 
Mean Temperature Differences (°C) and 
Sap memo candid (DeVid tO 1S is spe lsie%s eiaisie\ sss sss 65 
Cooling over the Yet Surface (°C) and 
Sapten Sanat OG FUE Vat LOM) shelslssiwicts <icis lcs s/ete a 1 


Stabilisation of the Temperature 
Gradient over the Cool Surface (°K 


m—~ 2) @enseeeeoeoeee2eaoeeeeegenvpeseeaeeeeesvesteeaeeeGBeeseeeees 79 


Wind and Diffusivity Parameters for 
FNPUG ee} the RPB Model eeeees2e28ees8se@2ene2eseseeevanaeeece828 84 


Experimental and Model Shape Factors ...... 87 


MoamiValicou ams eee Cen) 8 | set 93 
Z2 Sa Z, 
Upstream and Downstream Wind Profiles ..... She 


Vee 


J : at 
i ; 


, 


ni ; v 
Bes : Ry 


ae ‘ino 
f oa a a 
| i. Dy . wie Viv, iy are oo a ais 


» UT A ee Se 


To J aan) ; ; - yas 


whe 


= ow dd n abe meta covbreyawy bprne ria 
Bis. (2°) ADTIOTS =e smipaasnast pont . : +5 “ve. 
‘S sees eee a% et oan + pn Sebeae ttyeSn : Ae. oy 


ioe (9%) sonar ce¥ ote mae =e 
onl4 ep ey eoapenes, SRNR ETEO OSE ii SLORRS. - a 
a ix 
BOW TRTS Sqn? o4% Ye. Sat) sa ui eee e 
1) snotoevl too? <4 Gee 2 a0 cbaae a? 
bina be eae cabs tidedesskeuseean ae (*~@) : . 
— ; . ae 
36° Sot persed. eFiviews id - bis Debs Fae 
-¢ és a o's > &@ ees EG0% asa Qi z aed sagas ’ i : 
68 consen WLODOMN OGel® Leth buy iasnen ae 
‘ : me? eae a 4 saute’ eet 
. - > " 


svew Beto 029 Lah! @eesineeod ban igeitaqe :. 


TABLE Page 


LIST OF FIGURES 


FIGURE Page 
1-4 Two-Dimensional Internal Boundary 

LAY SU Mejaisie cere s co sm ose oe 0 0c sins 00 6 5 ssie sc e's oe 5 
3-1 Shwetz and Philip Solutions for IBL 

BKEISCNAe ous Son nic otioD OOCID A> DOO DOO Oo 6 domo ce 31 
el, SEVUCTUREROf VenperaturSe Sensors (0.6: s/s6- « 44 
6-2 SEVEKORE LO ie SG OOM A PASS 6 Ais 6 re Chou Pacacer poi Oia 45 
6-5 Pattern Of DIfLerences ..ss ccc scvcevcenecs 47] 
6-4 THERE XPeSEENe Med lr ote Gels o's ale alsicls oo 4 6 u)s,5's sale eas 
6-5 er Sic ISU eae I Any OU 4 its letelele ls efelalalstelelels ois ieisis sets 41s 60 
6-6 Photograph Showing Layout ...scccawocccccces 61 
Be Distribution of Temperature 

Dat ESRCICCSH iets cic cml ow ons soe wieisie cc us 0 0610's 0's 74 
<7 Cooling over the Wet Surface 2... .202cccecas 76 
hs Cooling over the Wet Surface under 


Stable andetinstable Stratrvrrcation sacs <c< aA 


7-4 Stabilisation of the Temperature 

Gradient over the Cool Surface .cccceceeene 81 
a Experimental and Model Shape Factor ....... 88 
J=6 Model Modification versus U, /K, --eesceeceee 90 
hoa | Model and Experimental Mean 

Stabilisation ee0ee@eesuveeseseeee2asegg#eseesteeeeeees2eeteesee2eseee?2e0 94 


saopone v4 re 
re 


1's 
a y 


fap #,% 


| boa 
rcahamoe’ Lantecat segodtnenko-ser 


 eacdb'vh Sesee dattbaasa tapbua ass enn woes 

1% dot eseotéaiet? ukelay viavde her 

p wele« de det ddechaldcd power ones deve e CoPaee ey ea 
move eset og . | ee 
of «Vee beeeun Jce YAU eee TOuneE 


tee svaeen ouee FaoAgia3. 2G 16 CTAIING - 


wee tae o@ & © G ed : ietineas i tears “ot? 
Yi 
4 
essere wee “toon @pease ova 7 POF Ne >nenutsaat 
soGqur pmbwode aqosposede 
ae we © ae@eweeoute ewe POIs ‘Cen fel Fe ai S2POTG 


aTrguscose? So wols ow daaeea 
** it oo) eee a ec oo 


ebenecewee O0ES502) 300 Sad aove gaiiocd 


sofa ecetaee gem 043 gava. pniloos a 
soe HOF PRVLTLIAIIE SLeaIsHl) Hus. Hidaze |, 


oypeessgee? Gis te goignehiidare 
(A -oessesveiee SOSRaME BGG? Gh? Tepe tnolhs38 


«08 dotnet ena aber dna sapehih doth ve 
ge!) se eden Ta aint nodsmatil bon Scant T 


bt? Te ae 


af fuih aor — 


FIGURE Page 


fale | Average Deceleration over the Cool 
Surface eeeaenreeseeveeueesesseeeoeeseespteaeaeeeeeeeee#sess8ts 8 @ @ 101 


<i 


i” y n : 
: 


| € fF * Merve “ees niece ae iy | - r : ; 


= 
_ 
— 
os 

a 

- 


CHAPTER 1 
INTRODUCTION 


The wind transport of any entity, such as heat, 
momentun, or moisture, which results from surface 
inhomogeneity on the nmicrometeorological scale is termed 
"local advection’. 

Although the phenomenon of local advection has 
been experienced by man throughout history, it has received 
considerable attention since about 1950 for several reasons. 

It has been realised that local advection is an 
important aspect of agricultural sScience.- For example, 
because the microclimate ree a finite irrigated area of 
land is affected by the advection of drier air from upstrean 
unirrigated land; consideration’ "o£f° local advection is 
necessary for the prediction of water needs for irrigation. 

Local advection over smaller distances results in 
interaction between experimental plots, as pointed out by 
Rider, Philip, and Bradley (1963) who stated that 

Frequently the dimensions of plots, subjected to 
different treatments, may be of the order of the 
range in which advective effects are demonstrably 
imoortant. ft 1S (clear that hetereogeneity 
between Such plots due to differential water 
treatments, differential aerodynamic roughness or 


differential radiative balance, may induce 
Significant unwanted interactions between plots. 


- 


It is clear that the value and applicability of 
plot experiments to large-scale field experience 
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can, in some circumstances, be very greatly 
enhanced if the presence of advective effects is 
recognised and appropriate allowances made. 
The development and testing of flux and profile 
laws for the atmospheric boundary-layer lead to increased 
consideration of the fact that airflow above the ground is 


very rarely in equilibrium with the surface immediately 


pelow ) it. Surface inhomogeneity in roughness, elevation, 


temperature, wetness, vegetative cover, soil type, etc., ~ 


prevents equilibrium because the influence of any 
disturbance is advected downstream more rapidly than it is 
propogated vertically, and only in airflow over an extensive 
horizontally-homogeneous surface will the boundary-layer 
characteristics be consistent with the surface boundary 
Eee One imposed. Conseguently, theoretical and 
experimental meteorologists began to study disturbed 
boundary layers, concentrating, in “the anterests of 
Simplicity, on the two-dimensional internal boundary layer 
(IBL), which develops within the boundary layer when flow 
over a uniform surface encounters a step change to new 
uniform, conditions, (see Fig. 1-1.) . Although attention 
has been concentrated on the influence of the surface on 
turbulence statistics (mean windspeed, turbulent kinetic 
energy, Shear stress, mean temperature) information gained 
Chee chem veftcc’ Pony “alerlow Of “a, Suridce discontinuaty 2s 
information gained on local advection. 

One of the most extensively investigated aspects 


of airflow over surface inhomogeneity has been the IBL 
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Figure 1-1. Two-dimensional internal boundary layer. 


& (x) height of IBL 
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created in flow over a step change in surface roughness. 
This is probably a consequence of the need to develop 
criteria for the exposure of anemometers, in order to 
ascertain that windspeed measurements to be related to the 
underlying surface are taken within the layer of 
re-established equilibriun downwind of the nearest 
discontinuity. 

Evaporation from a finite wetted area is another 
surface discontinuity, or local advection problem that has 
been widely studied. This problem is very similar to the 
problem of a surface discontinuity in temperature, if heat 
is treated as a passive atmospheric property. In nature, 
both discontinuities usually occur together, because an 
alteration in surface moisture conditions will alter the 
partition of available energy into sensible-.and latent heat 
fluxes, and thus create a auscontinurey in surface 
temperature. The moisture and temperature fields may be 
considered to form a field of ‘equivalent temperature', 
defined by ( L ) 

Cie eel TN oes (PS 


and under certain circumstances the finite wetted area 


advection problem may be fully formulated in terms of Can 
(Symbols used in this report are defined in an appendix. 
They are only re-defined in the text in cases where this is 
judged appropriate.) 

Sutton (1953) showed that the local evaporation 


rate decreases with distance from the leading edge, because 


of enrichment of the air with vapour. This indicates that 
surface temperature will not be constant downstream of the 
leading edge, so that horizontal homogeneity cannot be 
assumed for either temperature or moisture fields. In 
addition, because the temperature change may Significantly 
affect the longwave radiation ea Té, and moisture may 
affect both radiation and soil heat flux, all components of 
the surface energy balance equation may be changed by the 
discontinuity. 

In the experiment to be reported in this thesis, 
hereafter referred to as the ‘Ellerslie experiment', the 
primary aim was to study modification of the temperature 
poolLitiesataiter sa tistep change in surface temperature. 
However, as the surface temperature change was obtained by 
itrigating a patch of a large, reasonably. uniform grass 
field, the Boa set up was the complex moisture-heat 


interaction previously discussed. 
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CHAP OE Raz 
THEORETICAL BACKGROUND FOR AN INTERNAL BOUNDARY-LAYER 


2.1 GOVERNING EQUATIONS 
The basic equations governing Stratified airflow 
on the micrometeorological Scale are : 
The equation of continuity 
dA A Vipy a) (oes 
ot ~ 
The momentum equation 
dy = —1Vpb wae + yp Vy » Syne all (ee) 
dt p Je cares fe: 


The first law of thermodynamics, or ‘energy equation' 


+ 
pa rE bV-¥ =kV°T + p? = pao (ues 
The equation of state 


p a jae : (2.1.4) 


Here VeT, P,and p are instantaneous values of the velocity 

vector, temperature, density, and pressure, 0) is the 

instantaneous dissipation of kinetic energy into heat, and 
R is the rate of heat addition by radiation. 

A set of eguations describing the mean flow of a 


turbulent stratified fluid may be derived from the above by 
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rewriting the instantaneous equations using the Boussinesgq 
approximation, which assumes that density fluctuations are 
negligible unless multiplied by gravity, then splitting the 
flow into mean and fluctuating components using the Reynolds 
averaging convention. This procedure is demonstrated by 
Plate (1971). 
The resulting equations for the mean flow are, 

using the summation convention ; 
The equation of continuity 

Oui =a) (etl) 

Oncy 


The momentum equation 


du: | db al 1 du: 
LEST a eee OLE J Cae A bad (Zio) 
dt Peel deck: ag ee 


The thermodynamic eguation 


Ongar cme eta (2.1.7) 
At 3%; 


where a prime denotes a fluctuating quantity and all other 
quantities are now mean values. T, and fe are values for a 
hypothetical reference state with adiabatic lapse rate and 
hydrostatic equilibriun. 

In order to obtain these results viscous friction, 
molecular heat transfer and radiative flux divergence have 
been neglected. 

Equations (2.1.5)-(2.1.7), augmented by such 
additional conservation equations aS are necessary to 


describe a particular problem, are the basis for the 
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mathematical formulation of the internal boundary-layer 
problem. They. apply.- to Small-scale motion in the 
atmospheric boundary layer regardless of the presence of 
surface inhomogeneity, and are therefore valid inside, 


outside, and on the boundary of any IBL. 


2.2 Elimination of Velocity as a Variable 

In the analytical models to be discussed in 
Chapter 3, the number of equations is reduced and the 
problem simplified by the use of pre-specified velocity 
profiles, which are assumed to describe the flow both inside 
and outside the IBL. Through this approximation the 
variable W, becomes Simply a parameter in the IBL equations, 
and the momentum and continuity equations are no longer 
required. Given a realistic approximation for the velocity 
profite, “the ‘approach “shoul dawben) satasfactorye Pior the 
description. sof seany *IBLeset up by a discontinuity which is 
Not #ikelyetto, affects thes vedocity profile -(e.g.,. ‘step 
introduction of a conpletely passive contaminant). However 
in other cases , modification of the velocity profile is one 
of the important effects of the discontinuity (e.g.-, change 
ia sroughness)~ andesage full, formulation jof thei aproblen, 


including the momentum equation, is desirable. 
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Za3 Power Law Windspeed and Diffusivity Profiles 


In the theoretical models of local advection which 
will be described in Chapter 3 the wind profile is commonly 


approximated by the power law 


ie u, (Z| es @ (Zee 98) 


K = Ki (2) : e320} 


The true wind profile has been found to have the 


curvature 
PORUME, we yi, SMES ae z? du 
3(ln z)* 32 oz? 
>0 for Stable stratification 
=0 for neutral stratification 
<0 for unstable stratification 


(see Sutton (1953)). 


The restrictions placed on m are then 


m2>0 stable 
m2=0 neutral 
m2<0 unstable stratification. 
In order that u increase with height, m>0 is required. It 


can thus be seen that the power law fit to unstable and 
neutral wind profiles is not very accurate. 

The vertical fluxes in the lower atmosphere are 
measured to be approximately constant under the following 


conditions: 
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(1) Steady flow 
(2) Horizontal uniformity in all properties 
(3) Zero radiative flux divergence. 
Under such conditions a relation known as "Schmidt's 
conjugate power law' holds, i.e., 

n=1-n 
where n is the exponent in the power law for momentum 
diffusivity. If this law holds, measurement of the wind 
profile giveS an approximation to the momentum diffusivity 
profile, and through assumed relations between K, , K; , and 
Ke , an approximation to the other diffusivities. 

In the models to be described in this thesis, 
interest is in the wind profiles in the lowest 2 m of the 
atmosphere, where stability effects are expected to be 
unimportant, and the wind profiles should show the shape 
characteristic of neutral Stability regardless of 
Stratification. To represent a neutral profile we require 

ee oe 7 PS ay mamta lal zi aiarz 
From this formula, nm is undefined at z=0 or z=1, has the 
value 0.36 at z=10, and decreases with increasing height 
thereafter. 

The riiee of m for neutral conditions is commonly 


chosen as + , a figure originally based on wind tunnel 


5 
experiments. However, this value is most successful in a 
range of heights of 2 to 4m, whereas in the Ellerslie 


experiment m was chosen to best represent the wind profile 


below 2 m, and values of around 0.26 were required. 
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Under stable stratification m is expected to 
increase towards unity, while for unstable stratification it 
Should decrease. However, as has been noted, very close to 
the ground the wind profile is always approximately neutral, 
and in the Ellerslie experiment values of m were always 
between 0.24 and 0.30. These values agree well with those 
reported by Morgan, Pruitt, and Lourence (1971). 

Use of the power law relations and Schmidt's 
conjugate power law iS most convenient in analytical 
formulations, but it must be remembered that the technigue 


is only an approximation. 
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CHAPTER 3 


ANALYTICAL INTERNAL BOUNDARY-LAYER MODELS 


aoa The Model of Sutton 


~<a a 


The steady-state two-dimensional diffusion of heat 
in incompressible turbulent flow may be described by the 
equation 

Ue cel TS err (K ar | (B21. 1) 
Jx dz : 
This equation is obtained from the thermodynamic equation 
for turbulent flow 224.17) i modelling the vertical heat 
flux using the mixing-length approach 
Bee ene (3.1.2) 
dz 
Molecular heat transfer is not included, and the divergence 
of the horizontal turbulent heat flux has been neglected. 

Equation (3.1.1) may be generalised to describe 

the two-dimensional diffusion of any pasSive property in 


incompressible flow, 


U 36 ae G a (ac) 


02 Ree NS OO Ns 
where @ is the concentration of the atmospheric property 


under consideration. 


Sutton (1953) gave a solution of (3.1.3) subject 
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to the boundary conditions: 
Lin OX, Z)=0- 0 ToL x20 
Z-70 
Lim OX 7Z) =O FOr EZ > 0 
x0 
Lin 06 (x,z)=0l Horrex>. 0 
Zo 
If @ is regarded as a temperature difference T-T, 
Sutton's solution gives the temperature in the IBL that is 
formed when neutrally-stratified air at temperature fT, 
encounterS a step jump in surface temperature to a new 
Surface temperature T,. 
Sutton used power-law expressions for wind and 
da£fusivity. He related the two profiles through an 
approximation for the Lagrangian autocorrelation coefficient 


to obtain 


aci-n) 
K = a (2) 2-n 
Z, 
wieresra WiIsved. Krowneat unceronNrosean a, Z; and kinematic 
viscosity; Z is the reference height. For neutral 


cond veions, 1n=0225% 


The solution for temperature within the IBL is 


ies z)-To =(T.- TL) ft aisintae HaGaen | 


Ynere- bers a Known LUNCELONZOf UU, 72, 70s yy and: ={, is the 


incomplete gamma function, 
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BUT TiSTeintetesting to note that this solution is 
VEGyY Samtlar Gn yLOr malo, thatwontasnedar by Philip (1959) : 
(see Section 3.3), for the same problem, probably differing 


only in the choice of power-law coefficients. 


<j 


a2 The Microclimate Modification Model of de 


ries 


De Vries (1959) presented a theory of the 
influence of sustained irrigation on the energy balance and 
climate near the ground. 


The variables considered were 


Ons diatislan, 

(i ask 
whence ui" DebersetowrErigaveds@land@values and = td *orefers “to 
dry land values. Guided by measurements, de Vries neglected 
time changes in 8 and q, which he states are small compared 
toathesspatrval'changess “The "diffusion “equation, 9(35153) , Ys 
the basis of the analysis, and all variables are considered 
to be averages over several days. 


The wind and diffusivity profiles are represented 


by usu (i+ 2)" 
Ker 


and it is assumed that these profiles are unaltered as_ the 
air flows over the irrigated surface. 


The boundary conditions are: 
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(1) Lin q (x,z)=0 
Z—> od 
Lin 8 (x,zZ)=0 
Z—> Ad 
(ii) A Simplified time-averaged energy balance condition at 
the surface for temperature 
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(iii) A water balance equation at the surface 
ie (x, 24) ” ae 
ie Oe ie 


where I is the irrigation rate, assuming no runoff. 

The solution is complicated, and Rider, Philip, 
and Bradley (1963) stated that numerical solutions are very 
laborious to calculate. 

De Vries tested his model by setting up four 
climatic recording Stations: one over unirrigated land , the 
other three up to six kilometers downwind over irrigated 
pasture. Starting with the meteorological data for the dry 
land, averaged over periods of a few days or longer, the 
average temperature and moisture profiles over the irrigated 


area were calculated as functions of downwind distance. The 


predictions were proven to be reasonably accurate. The 
standard errors in the upstream-downstream differences 
were 0.3°C (in temperature) and OF 7°Civ(zn dewpoirnt)s 


Typical differences were: 
AT(Q cm) =5°C 
SEIS) cnyins icc 
AT, (0 cm) =6°C 


AT (125¢ca))—12¢ 
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The de Vries model also provides a means of 
eStimating the potential evaporation rate and thus the water 


needs of irrigated crops. 
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Rider, Philip, and Bradley 


Philip (1959) presented a theory of local 
agvectiom. In which §(3%1.3)s lustesolived for nconcentrationy 
flux, and radiation boundary conditions. 

Diegpowersvaw cquations "(2.3.1)" and (2.3.2) are 
used for windspeed and diffusivity!. Through this approach 
the variable u becomes simply a pre-specified parameter in 
the IBL equations and the momentum equation is no longer 
required. The theory is developed for any m and n, but in 
this (discusSionesi ti? owilLl be limited to the case where 
Schmidt's conjugate power-law holds, n=1-m. 

Philip's method is as follows. The transformation 


'- 2m 1+2m 
i = Uma EZ (3.3.1) 
(i+ 2m)* K, x 


reduces (3.1.3) and the power-law equations to the equation 


d* 8 dé [+m a) a (Bh3e 0} 
a nt ay dy ( | i rel m | 0 i 


It is straighforward to integrate (3.3.2) twice to obtain 


the general solution for the concentration of 9, 


teineriacts TPhilip. (1959) andwekRidergebhiliapy, and ibradiley 
(1963) used power-law forms without the z, which is included 
Tmet(Z soe icand. .(2Z.oee". The forms given in this report are 
dimensionally correct. 
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Pigetoes DtLerestings to note. thatthe parameters u, and kK, 
eer ecta ne solution onivs through: their ratio. u,/K, 7 which 
appears in fie 
Constants B and C are evaluated by using the boundary 
conditions, which nay be of three basic types: 
(i) Concentration Type 
Bes CS goo OEZ9(Z) %02-C- 850-0008 O(z)whenex=0,2>0¢. 
9 (7 ) =8:(x) when ype peas oc ueuwnhenex20pandaz=0% 
aot lLux’ Type 
Lin Le SAO Seer Holes 8=constant when x=0,z2>0. 


RL fied 
B= d(x) when x20 and 2>0, 


where 0 
Oe oar 
99 
py = Lim (-K38) 


740 (0) 
(iii) Radiation Type 
This type of boundary condition involves a_ linear 
combination of concentration and flux conditions at Yi =0 
(a= Come CxO, Z=0)). 

Now, for example, ait P one considered 
neutrally-stratified air in equilibrium flow over a surface 
at temperature To passing onto a surface with constant 
temperature T,, the relevant boundary conditions would be of 
the concentration type: 


ba Lely. 7 2 Os) Tat, when, x=0,270 
n> A 
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T=Eattor m= 0 i<e@.  T=T, when x20, z=0. 
Appisecatqon ofv these boundary conditions to (S-23.3))qives 
the solution for air temperature downstream from the 


discontinuity, 


[+m 
TiGcreap = (T, aie L(y ) a | (323.4) 
where (+P) represents the incomplete gamma function ratio 


I (4, P) a “ ay ci dy ff a etd, ° 


Philip gave the following expression for 


evaluation of the incomplete gamma function ratio: 


}+m 
ali C4 Ris aoe) (3.3.5) 


| 


Die nt ee pater 4 
|! Br HG) Cpe)” | 


where pane deine 
EOCENE 
and J’ is the gamma function. The series converges rapidly 
so that I can be computed quickly to any desired accuracy. 
Flux boundary conditions are treated by Philip as 
follows. 


Rewrite (3.1.3) uSing the flux-gradient relation 


to obtain 
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Differentiate with respect to z to obtain 
Ly eee i (4 a 
Ke nyse Nee ie oz. ° 
This equation is analogous to (3.1.3), with 
Omreplaced: by () 
n replaced by -m 
m replaced by -n (=m-1) 
u replaced by K-}! 
K replaced by u7~!, 


ies solved by the method ofmsolution of (3.1.3) to’ obtain 
B' ae a caon) e | Oca es 


where 
nies Ke eal a 
t (paar) due x 
The flux boundary conditions may then be applied. 

Rider, Philip, and Bradley (1963) considered the 
steady-state two-dimensional diffusion of heat and moisture 
in the lower atmosphere for the case of a step change in 
surface conditions. 


The equations are 


Ce ae gt (3.3.6) 
3 aioe rs (Ke L) 

Je pace avon 
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Ke = Kone Kemmtce(S:) che 


The upstream boundary condition is 
aC O.z.) and Ga-1e(0) 2) tony PooetOr Ze 0 (S232 3) 


The downstream surface boundary condition must be supplied 


by the equation for energy balance at the surface, 
& 
(Pm kauk. ~ ed 0 (0) + 6 = Ao LAE (5385) 


and an equation for the availability of water for 


evaporation at the surface, 
e,=¢,€T(x,0)) for x>0,z-0. (3.3.10) 


(3.3.10) assumes that water is freely available at the 
surface. 
By defining equivalent temperature 
Oce-ule Gexcbyipcrye 
and a newW variable 
OP Cae. 1 Oe (x, z) BC 0} 2) 


(293.6) ¥and (323. 7)¢canhbe wortten 


s@ Sane Me , 
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The upstream boundary condition becomes 
@,=0 for x=0,z20 ., 
Bquationse(3-329)sandu(ssse10)sforithe downstrean 


boundary condition are linearized to obtain a boundary 
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condition on the surface flux of @,. 

The problem is then soluble by the method of 
Philip, and a most important conclusion of Rider, Philip, 
and Bradley is that, to a high degree of accuracy, the 
solution is identical to that which would be obtained by 
replacing the downstream boundary condition, (S25 559) | ane 
(3.3.10), by the conditions of constant surface temperature 
fT) end constant absolute humidity e, =e. (Tp for, x>0. This 
is equivalent to the assumption that the quantities 

(am) Re +Ra - eg T*(x, 0) + G 
and A + L pu E 
are independent of x. Therefore it is concluded that the 
primary influence of the discontinuity in moisture and 
temperature conditions ieee change in the energy partition 
between sensible and latent heat flux. 

In eter tight ofMtenesesconclusions, Rider, 7 Philip, 
and Bradley Simplified their approach to the problem of 
modification of the mean flow by a step change in surface 
conditions. By defining 

[tFens Gers), -T(0,z) 
andu making suse “of the» fact’ that there is horizontal 
uniformity upstream from the discontinuity, the problem may 


be tCewrttten’ as the Solution of 
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subject to the concentration boundary conditions 
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% 
ar =e sud for x70,Z=0. 
The solution is obtained by specifying the general variable 
Garni. 3.3) equalgto.T “fuwith 
fh = U, wa Zan a 
(1+2m)* K, Be 
and using the boundary conditions on T* to evaluate the 


constants B and C. The result is the equation 
t+m 
Woz 1 Gun) = Gi-ss) Re ie Pomitioniy jisee (33 = 11) 


An analogous procedure leads to the equation for 


modification of the absolute humidity 


e(x,z)- e(0,z) = (.-6)| ACG EA CRO 


Rider, aglhi lip, ang Bradley gave a method of 
determining 
ATgatatggand WSS, Se. 
from Ke/K, »« me @(0,0) and T(0,0). 
With measurements of 
(i) upstream wind, temperature, and humidity profiles 
(ii) upstrean surface temperature and surface absolute 
humidity 
(211) athe, diffusivities,s whichignay bes «caiculated,..for 
example, from the relevant fluxes and gradients 
one can, under the assumptions of the theory, calculate 


downstream profiles of temperature and absolute humidity. 
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3.4 Limitations of the Rider, Philip, and Bradley Model 


There are many limitations to the model of Rider, 
Pot bese and@) bradley, (hereafter referred to as the 'RPB 
model'). 

As has been discussed, (Section Ps 8 ae the 
power-law profiles are generally a poor approximation to the 
real profiles. Furthermore, it is assumed that the wind and 
diffusivity profiles are unaltered by the change in the 
surface boundary conditions. While IBL situations in which 
this approximation could be accurate are conceivable, the 
case of a discontinuity in surface temperature and wetness 
foe noc such "a. -cace, “and mirc 25° nou likely that’ heat "can 
accurately be treated as a passive contaminant. 

That neglect of thermal stratification ‘iS “an 
important weakness in the model is confirmed by the 
numerical models of Taylor (see Section Wet) and Rao, 
Wyndaard, and “Cote™ "(sce Section 4.3). Taylor states that 
‘the importance of changes in thermal stability in the 
leading-edge problem should not be underestimated'. 

Downstream of any discontinuity there will exist 
hOCazZontaus gGradvents, —So, that “Che, fluxes) will) not) be 
constant with height and Schmidt's conjugate power-law 
connecting windspeed and diffusivity will no longer be 
valid. It has been pointed out that the RPB model is not 
restricted to the case where Schmidt's conjugate power-law 
TS. Vaadid. However, in all applications of the model 


discussed in this report the conjugate power-law has been 
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used, and the weaknesS in this approximation must be 
remembered. 

A further criticism of the RPB model concerns’ the 
assumption, embodied in boundary condition (3.3.10), that 
water is freely available at the downstream surface and that 
the absolute humidity at the surface is the saturation 
value. The numerical model of Rao, Wyngaard, and Cote 
showed that the assumption of saturation at the surface 
overpredicts the modification in the IBL. 

Fanallys 2 sshouldy, béeysanoted »ethat~athes user.of 
diffusivities is a device by which flux terms (covariances) 
arising in the boundary-layer equations (see Section 2.1) 
are related to vertical gradients. These relationships are 
known as 'K-theory! or the 'mixing-length model't and, since 
they define the diffusivities, their accuracy in expressing 
the fluxes .is-slimited; by, the), accuracy.; with, .which the 


diffusivities are specified. 
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220 The Model of Yih 


Yih (1952) found exact solutions to the general 
two-dimensional diffusion equation (3.1.3) uSing power-law 
expressions (2.1.3)-(2.1.4) for windspeed and diffusivity. 
The method used was dimensional analysis. 

When the diffusing entity is regarded as 
temperature and the power-law equations are specialised to 
the case n=i-m , Yih's solution is-the Same as that of 
Philip i(Section 3-23) for the flow of neutrally-stratified 
dumerhOn ea sUChacCe av. CONnstanu) tCemperavure {, Onto a surface 
at constant temperature T,. The two authors have solved the 
same set of equations subject to the same boundary 
conditions by completely different methods and obtained 
identical solutions, although this is not clear because of 
difference in nomenclature. 

Yih does not give a solution applicable to the 


case of non-neutral upstream flow. 
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SIis) The Shwetz Method 


Shwetz (1949) presented a method for approximating 
the solutions to several boundary-layer problems. His 
method has been summarised by Panchev, Donev, and Godev 
(1971), and may be applied to stationary turbulent diffusion 
-as follows. 

Multiply the two-dimensional diffusion equation, 


(3.1.3), by K(z) and introduce the new variables 


rf a ae (3.6.1) 
O(n) = u (a) K() 


to obtain 


Zz ° (3.6.2) 


Let the boundary conditions be 


(i) 8=1 at z=0,x>0 
(23) @=0 at x=0 for all z ('@9-less' upstream flow) 
(aa) Lin 8=0 

ZA 


corresponding to a step at x=0 from zero concentration of 4 
tou concentration “unity ) at the surface. These boundary 
conditions will necessitate the growth of an IBL, and the 
edge of this IBL is defined by 9(x,§ )=0. 

The value of 9 within the IBL is now given by 


(3.6.2) subject to the boundary conditions 
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Or ).= ; O(x,0) = (3.6.3) 


The Shwetz method is to find the solution with the 
help of successive approximations 
CeO TARO FTO. + 
where the individual contributions satisfy approximate forms 
Of (3-602) 5 
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(Shwetz gives no theoretical argument indicating that the 


UNE, a4. (8.605) 


above series approximation for 9 will be convergent. 
However, in his paper he shows that accurate solutions’ to 
many boundary-layer problems are obtained by this method 
even when the series iS terminated after the second term.) 

The first-order solution is obtained from (3.6.4) 
and boundary conditions (3.6.3) as 


RUE os 


The second-order solution is obtained using (3.6.4) 
“6, 
isa ik 


the expression 


Ge loe Ble [9994 oo 22 a 1 (6), » (3-6-5) 
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Application of the Peter ce conditions (3.6.3) gives 


(OVes : 0 } 


integrating twice, and adding 8, to give for 4 
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§ B 
sO, lay -| NES ( é 4) d 
a. Ae es dx 5 { y? peat P : (3.6.6) 


The double integrals in (3.6.5)-(3.6.6) may be Simplified by 
exchange of the order of integration. oa is evaluated by 
ae 


Stipulating that 


29 | = 0 (3.6.7) 
37] ula) 
and the equation for the growth of the IBL is 
o& ; SO aS ise (3.6.8) 
oS | he Sp ote ; 
1°) 


The second-order solution for the concentration in 


Be one vcueroundsto mie 
O©= |- n(%-/5z sles (3.6.9) 
wnere = = S$ yBly)ay 
T= (yy) ay 
Tz = §° 4? Bly) dy 
J, fs y PCy) dy. 


Equation (3.6.8) for the growth of the IBL may be rewritten 


ery bat 2 
a SB S ; (3.6.10) 
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ety. Comparison of Shwetz and Philip Solutions 


WOite uvand Kas power-law) functions off z, 


u =u (Z)" K=K(G) 


Assuming n=1-m (constant flux layer), the Shwetz parameters 


become Elis b Zi 


pp AIS ae 


Ziaie 20) 

*. U, bi 142m 
q = Ls (1+ 2m) - 
fie Meaui Died. 1Ee") 
‘ Zan Cie 3m) 

-2 
ah st U, b 1+3m 
‘ Z,°(1+ 3m) 4 
where L asia pees 


1+ 2m 


reese ae: 


The height of the IBL, given by (3.6.11), is 


“in i+2m 
vip ne (bs) = «| Function of (muK,z)]. (Sai) 
Now fore any SsethoLent, in,0K; 7 andeza,;4e(xpz) and Zs (x) as 
given by the Shwetz solution may be calculated. 
The solution given by Philip (see Section 3.3) for 
the case of 9-less upstream flow encountering a step jump at 


the surface from 9=0 to 9=1 is 


- 
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where 


1+#+2m 


5 : U, Zee 2 
K, C1+2m)* Xe ; 


The two solutions will now be compared in several different 
£20 esitudtions. It must be borne in mind that the height 
of the IBL will depend on its definition. For, othe vShwetz 
method, the definition is (3.6.7). For the Philip method, 
unless otherwise specified, IBL height will be defined to be 
the height at which the value of ® has decreased to 5 
percent of its surface value of 1.0. 

In the Ellerslie experiment the mean value of no 
OvVerrallyteials was 0.26 with a range of .0.25. to 0730.. “The 
Shwetz and Philip solutions will be compared for m=0.26, 
With a variety of values of the diffusivity. 

Under forced convection conditions 

K, =K,,=ku,y Zz where for a grassy surface with individual 
blades 5) cme in = Lengeh, u,=40 cme s—-!, This gives 
K, =K,, (30cm) =480 cm? st. 

For the Ellerslie data a considerable range of MK, 
values has been calculated, from 130 cm= s~* up to 830 cme 
oe Eee The corresponding values of u do not vary greatly, 
and their average is 180 cm s~! . 

Figure 3.1 compares the Shwetz and eke Md Bs bye: 


solutions. Two curves are given for the Philip model in the 


NewS 
7 
eh 
y Pp 
? 7 banacueaae eo" i ie va 
At: a at i 
riod oe. tenn co ae 
: _ t ii 
ag ma) 
' an, 
f 
4 
: : 
noi t ay hae Loos 
] 
. 
r 4 ~ 
eatin 
mI Me ay 
= 
4 i eis | j 
4 
.* 
g GPSp) of 
_ a 
. 1) SGF1  gbs 
. 
: . os pat . Pes , . : sa 
hy) 6 3 beplar ges hnoges 
est! Te Se eee a 
_ tr 7 . 54 : 
Be 3 =! at oa 
i an 


; 
os a rane 7 2 
ji G2 eS a , wid ane 


ch inhom gities as abs | nova by 


: ae & 4 
rh) eh ea ey 


3 


Height 
cm 
Shwetz 600 
200 
P 600 
Shwetz 400 
150 
P 400 
P 106 
100 
Shwetz 100 
P 100 
50 
200 400 600 800 1000 1200 1400 1600 1800 2000 x 
cm 


Figure 3-1. Shwetz and Philip Solutions for IBL Height. 
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case K, =100 cm2 s~!, the higher being the case in which IBL 
height is defined as the height at which ® drops to 0.1 
percent of its surface value’of 1.0. It can be seen that 
the height of the--IBL: depends: quite. strongly on the 
criterion used to define it. 

Although (Ser) 5.0) and (327.2) <apply to the 
artificial problem of '@-less' upstream flow encountering a 
step jump atythe surface froma! *O=0) ton G=1"; they werevfound 
to be useful in interpretation of the experimental 
temperature-profile modification data, because they give 


some idea of possible IBL height for given fetch. 
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CHAPTER 4 


NUMERTCAL INTERNAL BOUNDARY-LAYER MODELS 


4.1 The Mixing-Length Model of Taylor 

taylor, -(1970)) presented ay numerical model of 
airflow over a surface discontinuity in roughness, heat 
flux, or temperature. The model 2s restricted to 
two-dimensional steady flow. 

The following equations were used for the mean 
motion, and were assumed to be valid upstream and downstream 
from the discontinuity: 


(i) The equation of continuity (2.1.5) 


(ii) A horizontal momentum equation derived from (2.1.6) by 
neglecting the horizontal pressure gradient and 2 ul? 
terms 
OU dU Swe 
Nye «Wie PS eS Ye 


(iii) A thermodynamic equation derived from (2.1.7) by 


neglecting horizontal heat fluxes 


Si ay yA 
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The fluxes were modelled using the Businger-Dyer 


forms for the Monin-Obukhov Similarity functions; thus, for 


example, the —u’w’ or w2 term is specified by the gradient 
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ou 


eT given the appropriate similarity function. 


The formulation is, therefore 
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For neutral and unstable conditions, 


sont 
bn Sill ethan = tj ‘ (4.1.3) 


Pm 


while for stable conditions 
Zar a 
eee a aCe ear aee a) (4.1.4) 


It is hypothesised that buoyancy effects, e.g. in 
flow. /from.a cold to.a hot, surface, are all accounted for by 
a change in L, and hence a change in the intensity of 
Vertical, Mixing. »As Critical assumption 1s, that, the friction 
velocity, heat flux, velocity gradient, and temperature 
gradient are related by (4.1.1)- (4.1.4) not only in 


equilibrium flow above uniform surface conditions, but also 
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Within the IBL developing downstream from the surface 
discontinuity. 

No numerical results were obtained when the 
upstream flow was stable, because of problems with the 
finite-difference scheme. However satisfactory results were 
obtained for neutral and unstable upstream flow. 
Interesting predictions of the model were the following: 

(i) Velocity profiles were relatively, but not completely, 
insensitive to a step change in surface temperature. 
Largest changes in the velocity occurred at the outer edge 
of the IBL, with velocity increasing when flow encountered 
an increase in surface temperature. 

(ii) Treatment of heat as a paSSive property of the air 
gave important discrepancies (when compared with results 
with stability effects included) for downwind distances 
a t0 27 (100m if z,=1cm). At short downwind distances the 
errors were small, in agreement with the fact that at short 
fotcn vene (IBL- as, contined  “toy,ay shallow Layer in which 
Stability effects are expected to be unimportant. However, 
Taylor claims that stability effects become important as low 
as 25-50 cm and that in consequence simple heat-diffusion 
models will be inaccurate in cases where there is a 
significant change in the thermal conditions. 

The model was used to Simulate the experiment 
performed by Rider, Philip, and Bradley (1963), {see Section 
5.1), and the predicted profiles compared reasonably well 


with experiment. Taylor attributes much of the inaccuracy 
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in the numerical results to lack of data for input to the 
model. One of his most definite conclusions is that 'the 
importance of changes in thermal Stability in the 


leading-edge problem should not be underestimated. ! 


4.2  Weismants Extension of the Model of Taylor 

Taylor's numerical model has been extended by 
Weisman (1975) for application to the problem of the heat 
loss and evaporation from a body of water. 


A vapour continuity equation was added to Taylor's 


set of equations, 


aq OG peep anes He 
Hj he lelz PRPOZ 
where Err - pKe 3h 2 


The stability length was modified to include the buoyancy 
due to vapour stratification, becoming 
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The other modification chosen by Weisman was the 
use of a mixing-length that cannot increase without bound. 
This was found to be necessary to force the solution for a 
satisfactory distance downwind of the leading edge. He used 
the Blackadar mixing-length 
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which has the value } in the limit as z tends to infinity, 
aS opposed to the unbounded form k(ztz,) used by Taylor. 
Best results were obtained when 2 represented a height 
corresponding to the depth of the surface layer. 

Weisman concluded that the contribution of vapour 
buoyancy to sensible and latent heat fluxes and turbulent 
energy is significant, though less than the contribution of 


thermal buoyancy. 


4.3 The Numerical Model of Rao, Wyngaard, 
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Rao, Wyngaard, .,and*» Cote (1974) presented a 
Numerical model of airflow over a surface discontinuity. 
The model is restricted to steady two-dimensional flow. The 
equations used for the mean flow are the following, valid in 
all regions of the flow: 
(Aja newcduatizon of continuity (2. 125) 


(ii) A horizontal momentum equation derived from (2.1.6) by 


neglecting horizontal pressure gradient and 2 u* terms 
QU bhi Jukie. Seas 
ree a SZ 


(iii) A thermodynamic equation which neglects horizontal 


heat fluxes 
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(iv) A vapour transfer equation 
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where ? is the mean specific humidity and O, is the mean 
potential temperature. 

A higher-order closure technigue is used, giving a 
set of 16 partial differential equations for the model. The 
lower boundary conditions are based on the equilibriun 
flux-profile relations given by Businger et al. (1971) and 
Wyngaard et al. (197 4)%. 

The model was applied to flow from an extensive 
smooth dry area onto grassy wet terrain, thus simulating the 
experiment performed by Rider, Philip, and Bradley (1963), 
(see Section Sains Excellent agreement with the 
experimental data was obtained. Important predictions of 
the model are: 

(i) Most of the variation in the sensible and latent heat 
fluxes at the surface occurs within a few meters from the 
d2scontinuity. 

(ii) The assumption of Saturation at the surface 
overpredicts the modification occurring within the IBL. 

(iii) The neglect of the effects of thermal stratification 
becomes increasingly important at longer fetches from the 
discontinuity. Only for very: short (<5 mn) fetches can heat 
be treated aS a passive property of the airflow, because in 
this case the IBL is confined to the lowest layers of the 
atmosphere, where buoyancy effects are of lessened 
importance. 

At longer fetches, the neglect of thermal 


Shratarication. CeSULts in an eEroL. However, because the 
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temperature modification at given fetch increases with 
decreasing height, at low levels the error is ae small 
percentage ‘of .thev,actual* modification. Close to the top of 


the IBL the error becomes comparable with the modification. 


a4 The Numerical Model of Novikova 
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Panchev, Donev, and Godev (1971) reported briefly 
the work of several Russian meteorologists who have studied 
problems of surface inhomogeneity. 

Novikova (1969) gave an outline of a 
two-dimensional numerical model designed to study the 
effectsetof™ abediscontinuity in 'surfiace “'temperature ‘cand 
moisture on arbitrarily stratified upstream flow. The 
equations used were the continuity equation, the momentum 
equation, the thermodynamic equation, a moisture transfer 
equation, and the turbulent energy equation. FluxeS were 
modelled using the mixing-length formulation of Laichtman. 

Results were presented in graphical form, and are 
OGPELCeLe tuscan” thebycurrentrirstudy. Novikava gave no 


comparison with experimental data. 
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CHAPTER 5 


EXPERIMENTAL STUDIES OF TEMPERATURE-PROFILE MODIFICATION 


Diol The Experiment of Rider, Philip, and Bradley 


Rider, Philip, and Bradley (1963) presented the 
first substantial set of data on the local advection of heat 
and moisture. Under conditions of steady insolation, these 
authors studied temperature-, wind-, and humidity- profile 
modification downwind of a discontinuity from tarmac to 
freely- evaporating, closely-mown grass. The change in 


roughness at the discontinuity was estimated to be 


Zags AR hn Popes Tee wee wer 
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so that the change in roughness alone could be expected to 
produce important effects. However, the aim of the 
experiment was to study and compare with the Rider, Philip, 
and Bradley theory (see Section 3.3) the advection due to a 
discontinuity in moisture (and therefore temperature) 
conditions. 

Forty-three acceptable observations were obtained, 
each of 10 minutes duration, and profile masts were placed 
at 0,1,4, and 16 m from the discontinuity. In general, the 
temperature profile showed strongly unstable lapse rates 


over the tarmac with gradual development of an advective 
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inversion over the grass. 

The: Rider; Philip; cand yyBradley cmodely ofsSlocal 
advection was used to calculate the temperature- and 
humidity-profile nodification, the overall result of a large 


number of comparisons with observation being 
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Predictions were most accurate in cases of large AT, Ae, 
i.e. for low heights and large fetches. The measured 
temperature differences were used in the predictions for 
humidity modification, which were more accurate than the 
temperature predictions. Rider, Philip, and Bradley 
consider the discrepancy in temperature prediction to be due 
to the change in roughness which is not considered in the 
model. It is estimated that the diffusivity over the grass 
would be greater than that over the tarmac by the ratio 

K grass 2.7 

Ke Farmat 
and therefore that stronger mixing over the grass would lead 


to lower temperatures in the IBL than the model predicts. 


isaayo ? nt; 
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The increase in diffusivity through increased roughness 
would be partially countered by the increasing stability 
over the cooler grass, but the authors consider the effect 
of the roughness change on the diffusivity to be dominant in 
the air close to ground where stability effects are 
relatively weak. 

Priestley, in a discussion of the Rider, Philip, 
and Bradley paper (Rider, Philip, and Bradley (1965) ) does 
not agree that the poor temperature prediction is 
attributable solely to the roughness change. He points out 


that CoG Zak m>0 


) 


ee roan 
is a poor approximation in unstable stratification, which 
occurred in almost all the experiments. 

Rider, Philip, Gaad) Bradley concluded) that their 
model should give accurate predictions of advective effects 
caused in situations not involving a change in roughness. 
Their experimental data has been used for comparison with 
the predictions of most of the numerical models of local 


advection. 


 suiaseret hee 
‘hIaDEaTs yo nantes | 3 
Shh ae ae aalehe doitvnie society 
ad Shen TREE 908 PL kee TORS aide : 


an rs, a 
ais gives PEST) Re oe oe ie 


er bs . ‘ 
sai tiat gov bv nad. SR KOS: cece % ok oe: etme nie i: pe: 
woh ( (2akhy yeldiedar lini) qe tiokisy s540q poate 


i. 
= E pes eae Ee ai ana VW ie Pat lan Ai ae "eae 
vy ae 


, ~ 2 
- > . a ¥ ete 
s * 


tno Et#ciog- em Rates BAGS OHS 1G) A wae * aes 


. 


A 
4 


. 
~ 
A 


veut Meet “abnor ot) etl eee eee fess iz 
sy itenge e@2s ise Rama hoxiiamm 


- , 
? wa Pd Pele by aft * Keke a hy on ba a i 0 abi. ; : 7 - 
- * wes * ae r ' ? 7 > 5% ~ - éj ‘ i Mr a a s° +f ert re Sa ayip seu 
> 
<aeandDpuc wcntds @ pelewal- Joa -eonlteesean abel a 


> 


y pon tsaneee (Bebe Baap anet G4. ¢2 Ab (oteanl loyae ale 


ref 4s aivboee@ Leb bocene eeiose, Peed. FC toot soLbarg) 44 


43 


vations of Dyer and Crawford 


Dyer and Crawford (1965) reported data concerning 
airflow over dry unirrigated land onto a heavily-irrigated 
rye-grass field. Because the dry land received no rain all 
summer, while the rye- grass was supplied with 50 cm of 
water there was possibly a change in roughness at the 
boundary, but this is not stated. Measurements of 
temperature up to a height of 5 m were made as far 
downstream aS 200 m from the leading edge. 

Sensible-heat flux measurements at a height of 4 m 
and variable fetch indicated, on the assumption of zero heat 
flux at inversion height, that the heat flux at 4 m exceeded 
that at ground level by an amount which decreased with 


increasing fetch. In a two-dimensional problen, 
(eke ou 
A(o) - A(z) = Cp \ pe ya dz 
(0) 


so that the non-constant heat flux indicates that horizontal 
homogeneity in temperature had not been attained over the 


irrigated land. 
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CHAPTER 6 


INSTRUMENTATION, SITE, AND PROCEDURE OF THE ELLERSLIE 


EXPERIMENT 


6.1 The Temperature Sensors 

Wesely, Thurtell, and Tanner (1969) have described 
a fast-response unshielded resistance thermometer which is 
Subject to very small radiation error. The thermometers 
used in this experiment differ from the above only in 
SLtEuceuralve-detaris, and their structure is shown in Fig. 


pais 


solder 


painted wire frame 


Figure: 62156 Structure of ,Tenperature Sensors: 


45 


The frame was constructed of 1.6 mm diameter 
tinned copper wire, painted to form an insulated nount for 
the fine” wire,” deaving “a°”™Smnall*-spot’’ on * each’ side’ “for 
electrical contact. White paint was used in order that the 
equilibrium temperature of the wire would be as close as 
possible to that of the air. 8 to 10 cm of hey diameter 
platinum-plated tungsten wire were then wound on to. the 
frame, and the joints soldered. The resistance of the 
sensing element thus constructed was about 1500 ohms. The 


current supply circuitry was as shown in Fig. 6-2. 


to recorder 


potentiometer 


6.00 V 


Figure 6-2. Sensor Power-Supply. 


The current drawn by the element was of the order 


of 1 mA, and the voltage drop across it, V4, approximately 


1.5 Ve. A 5% change in the resistance of the fine wire would 
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Give a change in V,,+ of approximately 4 nv. 

Initially each sensing element had an individual 
6 V power-supply, but problems were encountered with battery 
Stability, and in operation all sensors were supplied by one 
regulated power-supply. 

Wesely, Thurtell, and Tanner conservatively 
calculated that solar heating could cause an air-wire 
cenperacure wditterence of 0. 15°Cjanwstibi~air, 0.099°C in a5 
Mm S-1 wind. Provided all thermometers are equally exposed 
to radiation, solar radiation should not seriously affect 
the measurement of temperature differences. Testing prior 
to the Ellerslie experiment indicated that radiative heating 


would not be a Significant problen. 


Gare The Temperature-Sensor Network 


it waS intended to Set up a network of 12 
resistance thermometers so that mean temperature differences 
within a spatial grid could be measured as accurately as 
possible. 

The signal from each of the sensors was carried by 
shielded cable to a Honeywell ‘electronik 15' multichannel 
recorder. The recorder stepped through a cycle of twelve 5 
second periods, during each of which the output voltages of 
a pair of sensors were differenced and amplified. The 
differencing pattern is shown in Fig. 6-3. 

Just before the end of the 5-second period the 


amplified difference was recorded. A thirty-minute record 
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Upstream Mast Ist Downstream Mast 2nd Downstream Mast 


temperature-independent reference 


Figure 6-3. Pattern of Differences. 


produced by the _ system contained, for each of the twelve 
temperature-sensor pairs; 30 differences each separated by 
one minute. Because the printing system was mechanically 
damped it is difficult to say whether the point plotted 
represented an average over the 5 seconds, or a shorter 
period. In the following discussion and throughout this 
report, each printed point will be regarded as an 
instantaneous value. 

Consider an experiment during which N values are 
taken for each temperature difference. From the N values a 
sample mean can be established. On the assumption that the 
population of differences is approximately normal (which is 
shown to be the case in Section 7.1) and that sampling was 


random, a Student's t-test can be used to estimate the 


} 


accuracy of the sample mean. Values of the standard 


ie “a na ‘ ar i ia at 
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deviation of temperature difference, as measured in the 
Ellerslie experiment , varied according to the spacing of 
sensors and the stability, the range of values being 0.04 to 
0.74°C. The t-test indicates that with a standard deviation 
of 0.3°C the chance of a 0.1°C or greater error in the mean 
temperature difference is 15% for N=20, but approximately 0% 
LOETN=120-Tetherchance oftanterron.oG80505°C ore greater = %is 
50% for N=20, 7% for n=120. The t-test indicates that if 
the desired accuracy in the temperature differences is 
better than 0.05°C then the value of N should be somewhat 
greater than 20. 

A simple criterion governing choice of the 
Sampling rate for mean profile measurements has been given 
by akarnaLs (197 5)its 


otx) \* 


> 
N > (—= 


where g{x) is the standard deviation of x, and W\xteis the 
desired accuracy. Unfortunately this criterion was not 
considered prior to execution of the Ellerslie experiment, 
in» which the ‘combination’ of ‘the need to ‘use ’a*single 
recorder for “many differences, the need to cover many 
positions with sensors, and the limitation on sampling time 
imposed by the need for stationary conditions restricted 
values of N to at most 40 (sunny conditions) and more often 
20, values which were erroneously thought to be 
satisfactory. In some cases N was further reduced by the 


rejection of data to ensure that averaging was performed 
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over stationary conditions. Confidence in the mean 
temperature differences could - have been improved by 
obtaining a greater numbersy,oL “data points) yany ‘each 


experiment. 


CALIBRATION 

Temperatures for the calibration procedure were 
measured uSing "ZEAL" mercury thermometers certified by the 
Atmospheric Environment Service and accurate to within 
LOO E. 

Calibration began by placing all 12 thermometers 
in a draughtless environment on a frame enclosed within a 
thick polystyrene box, the aim being to hold all sensors at 
an identical temperature. It. was established that about 20 
Minutes after being closed, the temperature within the empty 
polystyrene box would become uniform to within the accuracy 
of the thermometers available. However, with the frame and 
sensors in the box, great care had to be taken that there 
was in fact temperature uniformity. It was found that air 
contained in plastic film containers, which formed the 
weather protection for sensor wiring, would equillibrate so 
Slowly as to give an anomalous local temperature up to half 
an hour after closing the box. Holes drilled in the plastic 
cured this. 

When uniform temperature had been established, and 
thus all temperature differences were zero, the recorded 


differences were adjusted to chosen zero positions using the 
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potentiometers seen in Fig. 6-2. 

Calibration required measurement of the recorded 
shift away from zero position when a pair of Sensors was 
held at different (known) temperatures. To perform this 
measurement , use was made of a second polystyrene box which 
was ventilated with a flow of cold air. Even with the 
airflow the temperature in this box was found to be uniforn 
POMewLChin’ other Frlimitst fof taccuracyisfenthe calibrating 
thermometers. One by one each sensor was separated from the 
other 11 and the response (recorder divisions / °C ) was 
measured. The testS were repeated as often as time allowed 
in order to try to arrive at an accurate mean calibration. 
Initially it was not assumed that the calibration factor for 
all sensor-pair differences would be the same, but after 
many trials, the variation between repeated trials on any 
ohne pair waS Similar to that between pairs, so that a single 


number was eventually chosen for all pairs. 


BEHAVIOUR OF THE MEAN TEMPERATURE SENSOR NETWORK 

The resulting system behaved unsatisfactorily in 
only one respect. Trouble was experienced intermittently, 
with no channels being completely exempt, when the recorded 
number for zero temperature difference would jump or drift 
Tapialy + this phenomenon 1s ~hereatter referred to asia 
"Zero shift". This meant that the signal voltages arriving 


at the recorder were changing by some mechanism other than 


change in the resistance of the temperature sensors, because 
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any problem arising in the recorder would be expected to 
have Similar effects on all channels (this was not observed 
Bog Degathesrcase)~ Git was established thaththe zeroushifits 
were not caused by wetness of the various connections, but 
nonetheless equipment was stored inside to avoid wetting 
connections. With no explanation for this inconvenient 
behaviour being found, it was decided that the problen 
probably arose from changes in the resistance of some of the 
many joints and terminals in the system. Because the 
temperature sensor output voltages being measured were very 
Small, ~5 mV, small extraneous resistance changes would have 
a large effect. 

The problem was dealt with by frequent testing of 
the recorder zero positions, usually both before and after 
each experiment. This necessitated movement of all the 
Masts and cables, and transfer of sensors into a small 
volume in which mean temperature differences were assumed to 
be zero. Although this was very inconvenient, and may even 
have aggravated the problem, it was necessary in order that 
the data obtained Should not be misinterpreted. Tnestact, 
interpretation of the data was in many cases difficult and 


in several cases impossible, because of the zero shifts. 
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6.3 The Sonic Anemometer-Thermometer 


During each experimental run, vertical velocity 
and temperature at 280 cm were sampled at 50 or 100 Hz by a 
sonic anemometer-thermometer, filtered at 20 Hz, digitized, 
and recorded on magnetic tape. This data provided a means 
of calculating the heat flux, of interest in its own right, 
and also as an indicator of how steady the meteorological 
conditions were for the experiment. 

The sonic anemometer~-thermometer was of the type 
described by Kaimal and Businger (1966). 


Kaimal and Businger show that 


! 
Qe, im Ko, Gye am 2 
! Hide 2amTKf.T d 
Q, +KQ@ = i 


where 1 indicates properties of the upward pulse and 2 
indicates properties of the downward pulse; f is frequency, 
@’ change in Aphase, bwe Wchangeseinswerticainvelocity,eT 

change in temperature, c the mean speed of sound, and d_ the 
separation of the microphone-receiver pairs. The ratio of 
the frequencies is given by K=f,/f, . 

Recording is performed so that 

gd, y Kg, = X Ne 

Q +K@ =6 | 
where n, is the change in ny, the number recorded 
representing vertical windspeed, and ny is the change in 


n_, the number recorded representing temperature. 
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The sonic-anemometer record was used to calculate 
the senSible- heat flux, the covariance between fluctuations 


in vertical windspeed and fluctuations in temperature: 


A= Gare, (EE p AB n; =) Consianl <n l< 


Calibration performed several times throughout the 


Series of experiments determined & and BP : 


6.4 Shear-Stress Meter and Humidiometer 

The shear-stress meter consisted of three 
polystyrene propellor anemometers mounted accurately on a 
mast so as to be mutually perpendicular. The output of the 
anemometers was filtered, digitized, and recorded on 
magnetic tape. 

The anemometers had a response time of 0.3 S, a 
threshhold of 0.3 m s~1!, and a linear response to windspeeds 
from 1.2 to 22 m s~1. Low magnitude and frequent reversals 
Of Sethe tivertical wind lead to arcalibration problem? for the 
vertical anemometer. Berntsen (1971) considered this 
problem , and found that the non-linear propellor response 


has a more noticeable effect on heat-flux cospectra than on 


the spectra of the vertical wind. On average he found 


| Standard deviation of Ww Weenie = Pees 


(standard deviation of 1 pees Stress meter 
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Guae eg sii a 


(heat Flux) shear stress meter + fine wire 


The response of the propellors to the horizontal 
wind should be satisfactory except for very high 
frequencies. 


In view of the findings of Berntsen, quoted values 


a 
of u,=(-w VW) are twice the measured values. 


The flux of latent heat was measured using an 


Electromagnetic Research Corporation Lyman Alpha 
humidiometer. This was calibrated after the Ellerslie 
experiment was performed, prior to data analysis. Latent 


heat fluxes quoted in Chapter 7 were obtained by calculating 
the covariance between signals from the sonic anemometer and 
the humidiometer, which were mounted at 7z=280 cm and 
Separated horizontally by 40 cn. The large separation 
between humidiometer and anemometer probably lead to an 
underestimate of the latent heat flux, but this is not a 
critical point because the flux has not been used in any 


calculations. 
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6.5 Retrieval of Data from Magnetic Tape 

The magnetic tape recording system could record a 
haximum, OF 97S gchannels;. of ».which, «7 could be filtered at 
optional frequency before digitizing. The sampling rate was 
most often 50 or 100 Hz on each channel. 

The data were retrieved by a computer program 
which would calculate the means, standard deviations, and 
chosen covariances over any time period chosen. Covariances 
were Calculated using the formula 

ag - By Ry 
where a prime denotes a fluctuating guantity. 

Shear stress, the covariance between horizontal 

windspeed in the mean direction and the vertical windspeed, 


was calculated using the fornula 


vau Senta sigio Fasciatome Cos One nas Sta Wee 
where tanG= w/v , V is the horizontal windspeed in the 
mean direction, and OQ is the direction of the mean wind. 
Use of double precision ensured that variables did 
not overflow even when sums of products for covariances were 
accumulated for one hour of data. The program was tested on 


artacicraisdata. 
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6-6 Experimental Site 

The experiment was performed during late July and 
August of 1975 in a fairly level field at the University 
Farm, Ellerslie, Alberta. Pig 6-4 shows a plan of the 
field site. 

It can be Seen that the extent of horizontal 
ned ores upstream of the irrigated plot was at most about 
100 me. Had wind direction during August been predominantly 
Nu, as is usual, the site would have allowed the 
considerably longer uniform upwind distance planned for. 

The grass Was mown before experimentation began, 
in late July , so that roughness was uniform within the 
boundaries shown in Fig. 6-4, and of magnitude 
approximately 1 to 2 cn. 

Portable aluminum irrigation piping was available 
in 6.1 m lengths, and water was supplied to the system from 
a water truck. Use of the portable system allowed the water 
to be placed downwind of the 6 m tower for a range of wind 
directions of about 902. Six  ytengths “of ~piping: ~plus, fa 
perforated plastic hose forming the leading edge gave an 
irrigated area of about 20 m by 20 m, and the water truck 
generally placed about 1200 gallons of water on this area 


prior to an experiment. 
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Figure 6-4. The Experimental Site. 
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6.7 Instrument Layout 

The instruments were arranged to provide as much 
information as possible about the spatial changes in the air 
as it passed over the surface discontinuity and downwind 
over the irrigated plot. The fast-response instruments were 
placed at a height of 280 cm on the 6 m tower to monitor 
conditions in the oncoming flow. 

During several of the experiments, one of the net 
radiometers available was placed over the wet grass, while 
in ail experiments the net radiation over the dry grass was 
recorded. During each experiment two dewpoint sensors were 
placed on a small portable mast on the irrigated plot, and 
two were mounted upstream on the 6 m tower. 

The 6 m upStream tower held a set of propellor 
anemometers at heights of 30, 54, 114 ,233, and 470 cm, and 
temperature sensors at 30, 60, 120, and 240 cm. 

The aGTSt CORA masts were necessarily portable, 
because of variable wind direction. Two light 3 m towers 
were used, each holding four temperature sensors and three 
anemometers. The experiment was run for a period considered 
the minimum sufficient to give meaningful averages under the 
prevailing conditions ( 20 minutes if overcast, 40 minutes 
if sunny!) with the downstream masts at 5 and 15 m from the 


discontinuity. If the grass was still wet the masts were 


i As has been stated in Section 6.2, these averaging times 
did not yield a satisfactory number of data points in some 
experiments. 
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quickly moved to 10 and 20 m and the experiment continued. 
The guidelines provided by Rider, Philip, and 

Bradley (1963) were used to match the anemometer and 

temperature-sensor heights to the expected height of the 


IBL. Heights used were usually: 


FETCH (nm) HEIGHTS OF TEMPERATURE SENSORS (cm) 
0 30 60 120 240 
5 20 49 80 160 
10 20 40 80 160 
15 30 60 120 240 
20 30 60 120 240 


Figures 6-5 and 6-6 Show the layout of the instruments. 


6.8 Experimental Procedure 

In order that a good moisture discontinuity be 
obtained, the experiment could only be performed if there 
had been negligible rain for at least the previous 24 hours. 
If this condition was satisfied, and a favourable wind 
direction was forecast, the instruments were set up and 
tested early in the day. When conditions became 
satisfactory the irrigation piping was laid down and 
Watering began. It took about half an hour to discharge the 
1200 gallons, during which time the temperature sensors were 
all placed together and the recorder zeros established. The 


piping was then quickly removed (1/2 minute) and the 
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Figure 6-5. Instrument Layout. 
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Photograph showing Layout. 
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aqownstream masts placed on the plot ( 1 minute). Averaging 
then began, and was continued for 20 to 40 minutes after 
which time the downstream mastS were shifted and the 
performance repeated. The temperature sensors were then 
re-zeroed, and if conditions were stadt satisfactory 
re-irrigation began for another experiment. 

Throughout each day, wind profiles were measured 
continuously on the 6 a tower. This enabled the 
construction of a record of the upstream wind profile, and 
interpolations could be made for any desired time interval 
provided wind behaviour was reasonably steady. 

With the large number of instruments involved in 
the experiment, and the necessity of moving masts and 
cables, running the experiment required four people to be 


reasonably familiar with the procedure. 
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CHAPTER 7 


TEMPERATURE-—PROFILE MODIFICATION 


jie )ltemperature=Protaie Modification Data 

Six experiments, ranging in length from 12 to 39 
minutes and including a wide range of atmospheric stability 
conditions have been selected for detailed study. The 
remaining experiments were rejected because of unsuitable 
Wind direction, large time changes in atmospheric stability, 
or zero shifts in the mean- temperature network (see Section 
Gor ae 

Table 7-1 lists, for each experiment, the wind, 
cloud, energy balance, temperature, and stability conditions 
in the upstream flow. 

Table 7-2 gives the upstream and downstream mean 
temperature profiles observed in each experiment. 
Temperature differences between the masts are also given, 
but it should be noted that not all these differences were 
taken between sensors at equal heights, so that these 
differences were not all horizontal. Possible sources of 
inaccuracy in the temperature-difference data are, in order 
of importance: 

(i) Zero: shifts, (see, Section 6.2) wWhach lead to 


Systematic errors which were in many cases hard to 
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(a) Vertical 


240 

160 

120 

Height 80 
(cm) 60 
40 

30 


20 


EXPERIMENT NUMBER 1 


Differences 


0 


+0. 73.16) 
+0.42 (. 16) 
#0/otol(.91 0) 


To 


(cm) 
1000 


#123 5.120) 


+0.80(. 20) 


#0.51(.16) 


Tio 


(b) Differences between masts 


Table 7-2. 


T (1000, 30) -T(0, 30) = -0.29(.30) 


(2000 ,80)i-2 (0,30) iS = 0577135 6.8 3) 
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2000 


tile 22.20) 


n/a 


+0.53(.33) 


T20 


Mean Temperature Differences (°C) and Sample 
Standard Deviations. 


Day 


' T - (Vt 10.09! 
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EXPERIMENT NUMBER 2 


(a) Vertical Differences 


Fetch (cm) 
0 600 1500 
280 %. 
160 T; Tis 


12:0 +40.0/9:(223} 


Height 80 +0.63 (. 28) +0.44 (. 22) 
(cm) 60 +4+1.11(.48) 
40 #1.01(. 48) +0.70(.41) 


30 +2.41(.36) 


20 #1.42(. 44) +#1.05(.53) 


(b) Differences Between Masts 
T (600,160) -T (0,240) = +0.49 (.61) 


T (1500, 160) -T(0,240)= +0.44(.74) 


Table 7-2 (continued). Mean Temperature Differences (°C) 
and Sample Standard Deviations. 


(a) Vertical 


240 

160 

120 

Height 80 
(cm) 60 
40 

30 


20 


EXPERIMENT NUMBER 3 


Differences 
Fetch (cm) 

0 1100 
T, 

Ti 
+0.08 (.04) 

=) set. (ero) 
+0.10 (.07) 

=~ Om 0 (107) 
+O 2G) 

~0.52(.19) 


(>) Differences Between Masts 


Table 7-2 (continued). 
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2000 


Ta0 


+0.04(.06) 


-0.18 {. 14) 


—-0.40 (. 16) 


(1100, 160) -T(0, 240) = -0.06(.08) 


T (2000, 240)-T (0,240)= -0.03(.05) 


Mean Temperature Differences (°C) 


and Sample Standard Deviations. 


EXPERIMENT NUMBER 4 


(a) Vertical Differences 


— Renee See ee me ee ee 


Fetch (cm) 
0 1000 2000 
Za00 er. bey 
160 (ae 
P20 + OS ON SCS) 0.00 (.07) 
Height 80 —0.08(.05) 
(cm) 60 #402014 (. 13) zien Ore 11) 
40 =Onr C12) 
30, 7f0227-(238) ~0.28 (. 18) 
20 ~0.41(.17) 


{b) Differences Between Masts 


T (1000, 160) -T (0, 240)= +0.06 (.15) 


T (2000, 240) -T (0, 240)= -0.03(.09) 


Table 7-2 (continued). Mean Temperature Differences (°C) 
and Sample Standard Deviations. 
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EXPERIMENT NUMBER 5 


(a) Vertical Differences 


Fetch (cm) 
0 500 
240 -T. 
160 Ty 


120 -0.26(. 06) 


Height 80 -0.31(. 06) 
(cm) 60 =0.38 (240) 
40 -0.52(. 10) 


30 >> =0 589.5 13) 


20 -0.80(. 14) 


(b) Differences Between Masts 


T(500,160) -T(0,240) = -0.22(. 16) 


Table 7-2 (continued). Mean Temperature Differences (°C) 
and Sample Standard Deviations. 
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EXPERIMENT NUMBER 6 


(a) Vertical Differences 
Fetch (cm) 

0 1100 

2409 7 

160 a 


120 0-04 26:(406) 


Height 80 ~0-31 (#08) 
{cm) 6076-0 $33 (..09) 
40 -0.44(.08) 


30 -0.39(.16) 


20 -0.68(. 13) 


(b) Differences Between Masts 


T (1100, 160) -T(0,240)= -0.16(.12) 


Table 7-2 (continued). Mean Temperature Differences (°C) 
and Sampie Standard Deviations. 


70 


71 


detect, and could not be included in the uncertainty 
analysis. 

(i) An insufficient number of sample points in the 
Shorter experiments, leading to wide confidence 
Lintcs. 

The between-mast differences in experiments 1 and 
2 have not been used in all sections of the analysis, 
because they imply that warming occurred over the cool 
surface at several levels. A plausible explanation of this 
would be the occurrence of zero shifts in the sensor system, 
but it cannot be proven that such systematic error was 
responsible for the warming indicated. Such warming in 
experiments 1 and 2 would be inconsistent with the cooling 
over the wet surface which is expected and was observed in 
experiments 3 to 6. 

The cooling observed over the wet surface in 
experiments 3 to 6 is shown in Table 7-3, with an average 
cooling being given for each height on each mast. The 
contributions to each average from the respective 
experiments were weighted according to the length of the 
experiment. The uncertainties given are at the 95% 
confidence level, and were calculated by: 

{i) Calculating the standard deviation of the enlarged 
sample for each height 

(ii) Using the student's t-test with the above standard 
deviation. 


Formally, the student's t-test may only be used to obtain 
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confidence intervals for the sample mean when the sample is 
randomly selected from a normal distribution. In practise, 
however, it has been found that the confidence intervals 
defined by the t distribution are satisfactory even if the 
distribution of the variable of interest is only close to 
normal. 

Temperature differences were sampled randomly in 
the Ellerslie experiment. In order to test the normality of 
the distribution of the temperature differences a large 
Sample of 650 differences covering both unstable and stable 
conditions on the 26th of August was chosen. alge i Then 
shows the distribution of these differences on probability 
paper. It is evident that the distribution is close to 
normal, hence use of the t-test to derive confidence 
intervals for the temperature-difference data is justified. 
It has been noted in Section 6.2 that confidence in the mean 
temperature differences could have been improved by 
obtaining a greater number of sample points in each 
experiment. 

It can be seen from Table 7-3 that the amount of 
cooling does not differ significantly between the two masts 
at the 95% confidence level. If the cooling on both masts 
is combined to form a grand average, the resulting cooling 
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Temperoture 
Difference °C 


0.1 ] 2 5 10 20 30 40 50 60 70 80 90 95 98 «99 99.8 
Cumulative Frequency 


Figure 7-1. Distribution of Temperature Differences. 


1D 


240 cm -0.03+0.03 °C 
120 cm -0.06+0.03 

60 cn -0.23+0.03 

30 cn -0.40+0.06 ° 


The uncertainties above are again at the 95% confidence 
level, and are calculated by the method given earlier. Fig. 
7-2 shows this overall average cooling as well as that on 
the individual masts. 

If the downstream cooling data is broken into 


stable and unstable cases the results obtained are: 


Height (cm) Unstable Stable 

240 -0. 0340.03 °C n/a 

120 —-0.1140.04 -0.03+0.04 
60 -0.28+0.03 Oe te U 2106 
30 -0.57+0. 10 -0.284+0.08 


These values are plotted in Fig. 7-3. The cooling at 30 cna 
is significantly stronger under unstable conditions than 
under stable conditions, while at 60 cm and 120 cm the 
difference is not Significant at the 95% confidence level. 
Stronger cooling under unstable conditions is expected 
because in such cases the magnitude of the temperature step 
at the leading edge should be larger than under stable 


conditions. 
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Figure 7-2. Cooling over the Wet Surface. 
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Figure 7-3. Cooling over the Wet Surface under Stable and 
Unstable Stratification. 
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So far in the analysis of experimental results 
discussion has been concerned with the cooling of the air as 
rhe passes onto the wet surface, expressed as 
upstream-downstream temperature changes at each height 
considered. However the effect of the step change in 
surface temperature on the airfiow may also be considered 
from the point of view of the changes in the vertical 
temperature gradient. Such an approach does not require 
knowledge of the hetween-mast temperature differences, and 
thus experiments 1 and 2 may be included in the analysis. 

A temperature-gradient stabilisation parameter can 


be defined as follows: 


S(x, zm) F: iru ‘) [a2], 


Za ie Zi TED, a LZ 


It is expected that S will be positive regardless of the 
upstream stratification, and that it will decrease rapidly 
with height because a step decrease in surface temperature 
will have greatest effect on the temperature gradient close 
to the ground. Greatest values of S are likely to occur 
under conditions which maximise the change in surface 
temperature. 

Table 7-4 “dives a-the: (values -0of  S(x,7z)! fory abi 
experiments. In Vali but threes cases,).values “of S: are 
positive, the exceptions being only very small negative 
values. The expected decrease of S with increasing height 


is observed, and itS maximum values occur in experiment 2, 
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during which the net radiation, and therefore the step in 
Surface temperature, was larger than in any other 
experiment. 

It is also evident that the stabilisation of the 
temperature gradient relative to the upstream stratification 
increases with increasing fetch over the cool surface, but 
by an amount which is not significant at the 95% confidence 
level. Although it is to be expected that the downstream 
Seuratitication should be aincreaSingly stabilised with 
increaSing fetch until a new egquilibriun profile is 
established over the cool surface, the data indicates no 
Significant effect, and that averaging over both masts will 
now result in ae 1o0ss of information. 

PLC 7-4 shows the values of S plotted against 
height for each experiment. The mean value at each height 
is also shown. The error bars shown on the mean Curve give 
the uncertainties in the means at the 95% confidence level, 


calculated using the t-test. The mean values are: 


a1 02.CcH 050720207 °K nt 
113 0.1540. 06 
85 0.2840.11 
57 0.66+40.18 
42 1.33+0.32 


Because the temperature profiles were measured at 


varying fetches, the mean curve applies at the mean fetch, 
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Figure 7-4. Stabilisation of the Temperature Gradient over 
the Cool Surface. 
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Waech:) Was. 13°" A. The mean curve thus gives the average 
stabilisation expected at 13 m fetch over an irrigated grass 
surface of roughness 1 to 2 cm, with the averaging being 
applied over several experiments covering unstable and 
Staple stratification. | dans Section 7/-3 this experinental 
mean curve will be compared with that given by the Rider, 


Philip, and Bradley model of local advection. 


7.2 Calculation of Parameters Required by the Rider, 
Philip, and Bradley Model 
The model of Rider, Philip, and Bradley (1963) for 


temperature-profile modification is 


TT Cg) AAO, = Ab Cte Teas) 


[+ 2m 


where 
}- 2m 1+2m 


oe x 
Ge K, (1+ 2m)” XL : 


To compare the predictions of the model with the 
experimental data the parameters m, u,/K,, and AT, are 
reguired. The reference height z, was chosen as 30 cm, and 
T(0,z) was measured. 

The power law parameter m was chosen by plotting 
log u versus log z and taking the gradient between 30 and 
114 cm. 

The parameter K, was estimated from the heat flux 


(when available) uSing 
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ot 
where “yZ> 1s approximated by 


Se COMeO ae On) 


£0) 9130 

Under neutral conditions, assuming K,=K,, the 
formula K,;=kuyz gives a value for K,. This approximation 
has been used to obtain an estimate of K, for experiments 
during which the heat flux waS unavailable, but it must be 
remembered that values given by the shear stress meter for 
ux are not highly accurate (see Section 6.4). 

Table 7-5 lists the data available for diffusivity 
Calculations, the daftusivities, and the wind profile 
parameters. Accurate measurements of both temperature 
gradient and heat flux are required for an accurate measure 
of the diffusivity. Thus only in experiments 2 and 5 is the 
diffusivity regarded as accurate. In experiment 1 the flux 
is small and counter-gradient, while in experiment 6 the 


temperature gradient is not accurate. 
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7.3 Comparison of Observed Temperature-Profile Modification 


oe eee —S eee a ee et ne ee ae nee tree re ——— 


with the Rider, Philip, and Bradley Model 


ATo, the change in surface temperature, was not 
measured. In order to compare the model with the 
experimental results without the necessity of supplying an 
arbitrary value of AT,, the following shape factor was 


defined: 


Reha [T(x 2:)-T (x, z,) | brea [ T(0,z:) - T Co 2) lest 
[ese ec Z3) te Far Lege eG. Z:) Jexp 
where 2Z2>Z;>z,.- The model can be used to predict values of 


this shape factor. Define 
iam 
ec 2) = Custos & l+2m )) : (7.3.1) 


Then 
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The shape factor, model or experimental, can be interpreted 


physically by the following rearrangement: 
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is the temperature-gradient stabilisation in the lower layer 
with respect to the overall stabilisation through the deeper 
layer. 

For each experiment the F(x,z) were calculated 
uSing the appropriate values of u,,m, and K, (see Table7-6) 
in the formula given by Philip, (3.3.5).With the results of 
these calculations ripe is possible to compare the 
experimental and model shape factor values. Comparison by 
means of the shape factors solves the problem that neither 
the change in surface temperature nor (in experiments 1 and 
2) the horizontal changes in air temperature are available. 

Table 7-6 gives, for each experiment, the 
experimental shape factor and its uncertainty at the 95% 
confidence level, versus model prediction. Fig. 7-5 shows 
a plot of model and experimental shape factor values for all 
experiments. 

The values of K, used aS input to the model were 
those calculated from the heat fluxes, except for 
experiments 3 and 4, in which cases they were the values 
given by ku,yz. 

It can be seen that in all cases where the 
upstream temperature profile is stable (experiments 1, 5, 6) 
the model shape factor values exceed the experimental 
values, while the opposite is true when the upstream 
temperature profile is unstable (experiments 2, 3, 4). 
However, there are significant differences at the 95% 


confidence level only in experiments 1 and 2. In the case 
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Figure 7-5. Experimental and Model Shape Factor. 
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of experiment 1 the diffusivity used as input to the RPB 
model as ‘suspect. A larger) value “would seem more 
reasonable, and it is later shown that a larger value would 
give a better fit of the model and experimental shape factor 
values. 

It thus appears that the model overestimates the 
Shape factor in stable stratification and underestimates the 
Shape factor in unstable stratification , but the accuracy 
of the shape factor and diffusivity data is not sufficient 
to test this hypothesis stringently. 

The model values of the shape factor are dependent 
upen yithe “value Wot. (u,/K; - EPG). 7-6 gives a family of 
modification t.curves with different values of u, /K,, but all 
for x=10 m, m=0.26. The quantity plotted to represent model 
ROdntacationy ais -—E(e, Zz), detinedtby (7.5.1) and related to 
the modification because 

TRC a Coe RA ao) 
IE can be seen that as u, /K, 1S increased 
(i) the height of the IBL decreases and the difference 
between downstream and upstream temperature at any 
height is reduced. These effects are intuitively 
expected as windspeed increases OLea. fuse vi ty 
decreases. 
(a) the value of 
F(10, 60) - FC1o, 30) | 
F(10, 240) = F(10, 30) 


which is simply a model shape factor as previously 
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Figure 7-6. Model Modification versus u, /K,.« 
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defined, increases. 
It is evident that the larger the value of u,/K, , 
the larger will be the model shape factor values; i.e. the 
stronger will be the stabilisation in the lower layers with 


respect to that in the deeper layer. 


STABLE UPSTREAM 

The RPB model predicts excessive stabilisation in 
the lower layers, Rm>Rep - If a larger value of K, was used 
in the model, the shape factors would agree more closely, 
but upstream-downstream temperature differences and the 
height of the IBL would also be affected. In addition, one 
would expect the value of K, over the cooler surface to be 
smaller than that measured upstream, so that the use of a 


larger K, cannot be justified. 


UNSTABLE UPSTREAM 

The RPB model underpredicts the stabilisation in 
the lower layers, Ry<Rexp - There is some justification for 
the use of a smaller value of K,;, because the diffusivity 
over the cooler surface is expected to be lower than that 
upstream. However, a decrease of K, from 670 to 150 cm? s~! 
is necessary to match R, and R,, in the case of experiment 
2 This seems an exceSSive decrease, and would result in a 
very low height of the model IBL, as well as_ small 
upstream-downstream temperature differences. 


The experimental and model shape factors for each 
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height were averaged over all experiments, with the data 
from both masts being combined. The average values were 


found to be: 


Height (cm) Rexp Re 
40 0.3040.14 Q220 
60 0.53+0.11 On o2 
80 0.66+0.11 0.70 
120 0.8020.12 0.87 
These values are plotted in ios eG FS 7-5. 


Uncertainties are at the 95% confidence level. There is no 
Significant difference between experimental and model 


average Shape factor. 


7.4 Comparison of Model and Experimental Stabilisation 
In Section real a temperature-gradient 


stabilisation parameter has been defined and discussed. The 


RPB model can be used to predict values of this parameter: 
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where F is defined by (7.3.1). 
Because AT, was not measured, it was decided not 
to attempt experiment by case-by-case nodei-experimental 


comparison of this parameter, but rather to compare mean 
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curves. The model average curve for stabilisation versus 
height was calculated by calculating the values of 
F (X,Z2)-F(x,z,) for each experiment, uSing the appropriate 
NaLUCS Obs elle, Ky gue 2) and 9 iX, (see Section 7.2), then 
averaging over all experiments to obtain a mean value for 
each height. The values obtained are given in Table 7-7, 


and apply at the mean fetch, 13 on. 


F(z2) oa F(z.) 


TAPES say Se 


Zionteom) 


170 -9.05 
13 -0.09 
wehs, ! -0. 14 
57 -0.26 
42 = Onid 1 


F(z2) = F(z,) 


Z2 ro Za 


Table 7-7. Mean Values of 


Now if a value of AT. is chosen, mean model 
stabilisation may be compared with the experimental mean 
stabilisation. Fig. 7-7 compares model and experimental 
stabilisation for AT,= -2°C and AT, = -4°C. The shape of the 
model curve is quite close to that of the experimental 
curve. In view of the uncertainty in the experimental 


cutve, the model curve for AT,= -29°C does not differ 
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Figure 7-7. 


Model and Experimental Mean Stabilisation. 
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Significantly except at the lowest two levels. As has been 
discussed in Section 7.3, the model overpredicts the 
stabilisation at low levels with respect to the overall 
Stabilisation when the upstream flow is stable, and vice 
versa when the upstream flow is unstable. In the formation 
of a mean model curve, there has probably been some 
cancelling out of the over- and underprediction in 
individual cases, but the mean model curve nevertheless 
underpredicts the stabilisation in the lowest layers. 

The model curve may be shifted to larger or 
Smaller values of the stabilisation by choice of AT. In 
order that the ae curve should reproduce the greatest 
value of stabilisation observed at the lowest level, that of 
experiment 2 (net radiation 49 mW cm-2) the value AT, = -9°C 
is required. To reproduce the minimun lowest level 
Stabilisation, that of experiment 6 (net radiation 3 mW 
cm-2), the value AT,= -0.5°C is required. 

Thus the model mean curve indicates that 
throughout the six experiments the change in surface 
temperature at the leading edge varied between -0.5°C and 


-9°C, 
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CHAPTER 8 


MODIFICATION OF THE VELOCITY PROFILE AND THE RADIATIVE FLUX 


8.1 Experimental Wind Profiles 

The numerical model of Taylor (see Section 4.1) 
predicted that the velocity profile would be relatively 
insensitive to a step change in surface temperature. For 
example, for a very large step in temperature the velocity 


modification predicted is: 


These numbers ace taken from a graph given by Taylor (1970), 
and apply to the case 

Nios bt ocG 

x= 10 to 20 a 
given that u,=0.44 m s~t, the average value of u, observed 
in the Ellerslie experiment, and 2Z= 1 to 2 cm, the 
roughness length of the Ellerslie site. The maximum 
windspeed modification occurs in the range of heights 20 to 
40 cm, although this can only be estimated roughly from 
Taylor's graph. 


An attempt was made in the Ellerslie experiment to 
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test for velocity-profile modification and compare it with 
that predicted by the Taylor model. 

As mentioned in Section 6.8, wind profiles were 
measured continuously on the 6 m upstream tower, enabling 
construction of a record of the upstream wind profile. The 
anemometers used on the upstream tower were ‘'Rimco! 
metal-cup anemometers, calibrated to C.S.1-R.O. Standards. 

Downstream wind profiles were measured using 
'Rimco' mini-cup anemometers. The calibration of these 
anemometers was inaccurate below about 1.2 m s~!t , but only 
rarely were such speeds encountered. 

Possible sources of error in the wind profile data 


are the following: 


(1) Imprecise calibration of the downwind anemometers. 

ea) Use of different anemometers upstream and 
Poenet renee 

(111i) The interpolation required by the fact that upwind 
and downwind profiles were not measured 
Simultaneously. 


Table 8-1 gives the wind profile data, the 
upstream profiles being interpolated from the continuous 


record to give values at the times of the downwind profiles. 


(i) 7th August 
The trend in upstream windspeed, (a steady 
decrease), waS duplicated downstream. Both the 1410-1420 


and the 1424-1434 runs indicated that windspeed decreased 
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Significantly at all heights with increasing fetch over the 
cool surface. The average decrease over the 10 m_ between 
the two downstream masts was 0.3 m s~!?t. The average 
decrease over the 10 m between the upStream and the first 


downstream mast was approximately 0.14 m s~!. 


(ii) 14th August 

From 1326 to 1350 the upstream windspeed increased 
Slightly, by an average of 0.1 m s-t over all heights. The 
windspeed during the 1339-1349 downstream run Was an average 
of 0.14 m s7! lower than during the 1322-1336 run. Thus the 
upstream trend in windspeed was not duplicated downstrean. 
No explanation for this has been found. 

During the 1322-1336 run there was a Significant 
difference in windspeed between the 5 m and 15 m downstream 
masts at all heights, the speeds on the 15 m mast being 
lower by an average of 0.15 m s-!. The upstream windspeed 
was lower, though not significantly, than on the 5 m mast 
EOE Z245 cm, but was higher than on the 15 m mast for all 
levels. This suggests that the interpolation for upstrean 
windspeed may have given too low a value. 

During the 1339-1349 run windspeeds on the 5 o 
mast were significantiy lower than the interpolated values 
on the upstream mast, with the average value of the decrease 
being 0.25 m s~!. Windspeeds on the 15 m mast were lower 
thandson othetec5) va gaastirbys anitaverage .ofe0.1i ness), c20t 


significant considering the uncertainties. 
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(iii) 26th August 

During the 1623-1631 and 1634-1642 runs windspeed 
on the first downstream mast was significantly lower than on 
the upstream mast below about 50 cm. At greater heights the 
difference was not significant. Windspeeds on the two 
downstream masts did not differ significantly. 

During the 1647-1657 run the top anemometer on the 
first downstream mast indicated a significant increase in 
Windspeed in comparison with the interpolated upstream 
profile. At other levels, the differences were not 
SEGnNpe cant. 

During the 1659-1709 . run windspeeds on the 
downstream masts did not differ significantly. However the 
upstream (interpolated) windspeed was Significantly lower 
than downstrean for heights above about 60 cm. 

For the remainder of the afternoon of 26th August 
windspeed was variable, and the accuracy of the 
interpolations for upstream windspeed is dubious. In three 
of the five remaining runs windspeed downstream was higher 


than that upstream. 


Because the velocity-profile modification varied 
considerably from run to run, it was decided to average the 
deceleration at each height over all runs, putting together 
the data on the two downstream masts. The resulting average 


deceleration Curve 915. ‘shown an Fig. 8-1.) “fhe error bars 
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Height 
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Deceleration ms~! 


Figure 8-1. Average Deceleration over the Cool Surface. 
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given are at the 95% confidence level and were calculated 
using a paired t-test. Significant deceleration over the 
cool grass is indicated, for heights less than about 1 am. 
Because windspeed at ground (zero) cannot be modified, the 
maximum modification must have occurred between 35 cm and 
ground level. 

The magnitude of the observed average deceleration 
is inconsistent with the Taylor model, which indicates that 
a very large temperature discontinuity is necessary to give 
a Significant change in windspeed. The observed height of 
maximum deceleration does agree with the model estimate of 
205 LoAnlOnecmevLor se Letchest 1 o0fFG10 tor20um. Because of the 
possibility of systematic error in the wind profile data, no 
definite decision can be made as to the accuracy of the 
Taylor model. However, it can be tentatively concluded that 
the model underestinates velocity modification. 

A decreased horizontal windspeed over the cool 
grass implies by continuity the presence of a mean upward 
vertical velocity, or a mean difluence of the airflow around 
the irrigated area, or both of these possibilities. The 
existence of a mean upward vertical velocity would give 
large vertical fluxes. For exampie, consider the decrease 
in horizontal windspeed to be compensated for solely by a 
vertical component of windspeed over the wet grass. Using a 
value for the deceleration of -0.2 m.s-? in 10 m for all 
heights in the layer from 134 to 44 cm, the mean vertical 


welocity Sinéithe layer®isvoferhetorder of.ivcm sctc snehta 


“ew 
—_—" / 
sr) 
7 
an 
acs 


, s Ret ere & 
woLese 


ne me? 


103 


vertical velocity leads to a sensible heat Linx of 
approximately 300 mW cm-2 over the wet grass, and would also 
lead to a very large latent heat flux, depending on the 
absolute humidity above the grass. Hence the evaporation 


rate near the leading edge may be very large. 


Cer A Simple View of Wind-Profile Modification 


On the 7th August both wind profile experiments 
displayed significantly decreasing windspeed with increasing 
fetch over the cool surface. Particular confidence is held 
in these two experiments because of the consistency in the 
upstream and downstream windspeed trends. - In addition, 
Significant deceleration is indicated in the average over 
all experiments. That windspeed modification should be 
observed at all can be understood by the simple hypothesis 
that as the air passing over the cool surface becomes more 
stably stratified, mixing of momentum from higher levels 
becomes less efficient. However, the question is raised, 
'is it theoretically possible to have a difference in 
windspeed between the dry and the wet (cool) grass of the 
erder of 0t2antisathic®. 

An approach to answering this question is to 
choose a velocity-profile formula which incorporates 
stability effects and to see what velocity change is given 
by a change in stability. It must be remembered, however, 
that wind profile formulas are constructed using data 


gathered at sites where horizontally uniform and _ steady 
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conditions prevail, and that they are really only intended 
to be applied in such Situations. 

The change! in ‘stability (Ge...) change in LL, the 
Monin-Obukhov length) depends on the change in the heat 
flux. Unfortunately, no estimate of this can be made for 
the 7th August, when the heat flux upstream was small and 
counter-gradient. However an estimate of the change in 
stability can be made for the 14th August, on which day 
decreased Sneha over the cool grass was also observed 
although confidence in the accuracy is less. 

On the 14th’ August the temperature profile was 
unstable both npereoee and downstream. Businger (1973) gave 
the following expression for the wind profile in unstable 


stratification: 


Us {+x beoeiegy 


U= —& [he 22652 ae ees + Qarctan(x+Z)] 


‘ 
where aos CaGuas Vise e 2 Ne ° 
je 
A change in stability will lead to a change in u, 


assumed to be given approximately by 


it = NL ‘ 

Measurement of the upstream temperature profile 
and heat flux on 14th August allowed calculation of a heat 
diffusivity, K(30 cn) =670 cm2¥s-*.) Inyvorder «to» calculate 
the change in the heat flux (and hence in L) the diffusivity 
measured upstream was used with the downstream gradient on 


the 15 m fetch mast to give an estimate of the downstrean 


heat sctiuzxs Results were 
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For a height of 60 cm, the change in windspeed 
from upstream to downstream given by this simple method is 

u (downstream) -u(upstream)= -0.1 m s~!? 

Windspeed modification of this order of magnitude was 
observed on the 14th August. It thus seems possible that 
the considerable velocity-profile modification indicated in 
some of the Ellerslie experiments did occur, and that the 
Taylor model underestimates Such modification. 


8a Modification of the Net Radiative Flux 


In none of the experiments was a Significant 
difference shown in the mean net radiation between the wet 
and the dry surfaces, but in all cases the (insignificant) 
change was such that the net radiation was increased over 
the wet surface; for example 
(1) Experiment 2 50.2+0.8 mW cm-2 (wet) 

48.7+0.8 mW cm-2 (dry) 
(ii) Experiment 4 16.2+0.8 mW c,-2 (wet) 
14.8+0.8 mW cm-2 (dry 

However just prior to experiment 2 when a net 
radiometer was moved from dry to wet grass its output jumped 
by 4.441.5 mW cm-2, indicating a Significant increase in the 


net radiation over the wet grass. 
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The net radiation can be written 
R:, = (I-r) Rs + ike. tan caalie 
The materials comprising the earth's surface have 
emissivity values very close to 1.00 (see Gates (1965)). 
Typical values are: 
plants 0.97 
soils 0.95-0.98 
snow 0595-02995: 04. 
If a change in the emissivity was supposed to have caused 
the observed change in the net radiation between the wet and 
the dry ground, then the change in emissivity required is 
given by ee ee Bin : 
TT. 


Values of AR, and € using 1,=290°K are: 


ARn (mW cm-2) 


On the basis of these numbers, it seems unlikely that a 
change in the emissivity lead to the change in the net 
radiation. 

It is generally assumed that horizontal changes in 
R, and Rs, are negligible for a small irrigated area. 
However there was possibly a change in the reflectivity, fr. 
Kuhn and Suomi (1958) tabulate measured vaiues of cr over wet 
and dry . grass, but there is so much variability in the 


numbers given that no conclusion can be made. 


- = ny 
\ 4 joe Uta » 
‘ ' i a i ¥ : 
ges 7 CLO 
> ) eyes 
man 
rT i) ‘es a 
y . vey) 
1 7 ree : : ay a hain 
. ! . ae lee hi 
i} a ed d 
7 : ae 
es es ae 
@ 7 7 
a Cae We 
a eG 
* ms swell 4 
.' 
: 7 
: " Dn 
<a “s,s 
' > , 
. 
* 
en * . a ¢* 
. . « 
od = * 
er eas ea 
* * * 
. ey 
a ny 


A aphcpmtle+ttprgneeiiplalin release pet lag 
7 : Le 7 


— va! : : p. 
fa=% ? , ) 
or lee ee vey 
a } | 7 ihe 


¥ 
' i> P 


aia ol wok wat 


oan 


107 


The simplest explanation for the observed change 
in the net radiation is the change in surface temperature. 
A change in surface temperature gives the following change 


in the net radiation: 


INR = - eg ee AE 


Values of AR, for given AT, are ;: 


ARa (mW cm-2)) 0.5 

It was found that a change in surface temperature 
Of *¥about) *=9°C 7 Sis frequiredh® to match the RPB mean 
stabilisation to that observed during experiment 2 at 30 cm 
(see Section 7.4). As can be seen from the numbers above, 
such a change would give a change in the net radiation of 


approximately the size observed. 
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CHAPTER 9 


CONCLUSION 


oea Summary of Findings 

The local advection model of Philip (1959) has 
been compared with an extension of an approximate solution 
due to Shwetz (1949) for the artificial problem of '@-less' 
upstream flow encountering a Step change at the surface from 
8=0 to O=1. The two methods were found to’ give similar 
predictions for the growth of the IBL. The Shwetz solution 
leads to a Simple expression, Eq. (3.7.1), for the height 
OfVthe VIBLV@as-a~ functionsofefetchs 

From the six temperature-profile modification 
experiments analysed, mean results have been given for the 
cooling and the stabilisation of the temperature gradient 
over the wet grass. No significant difference was found 
between the cooling on the two downstream masts (separated 
by*10eny but ato=thee@lowest Picvel,” 30> cH; cooling was 


significantly stronger under unstable conditions than under 


stable conditions. Greatest stabilisation of the 
temperature-gradient occurred in the lowest layers 
considered. 


The Rider, Philip, and Bradley (1963) model of 


temperature-profile modification has been compared with the 
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experimental data of the Ellerslie experiment without the 
necessity ob estimation <of#eATs, ifthe) (change oinwsunface 
temperature, by definition of a shape factor (see Section 
Vis 34s Confidence intervals for the experimental data are 
Wide, and significant differences between model and 
experimental Shape factors were observed in only two 
experiments. A tentative conclusion is that the RPB_ model 
overestimates Stabilisation in the lowest layers with 
respect to that over a deeper layer in stable conditions, 
and vice versa in unstable conditions. The shape factors 
were averaged over all experiments, and there was found to 
no significant difference between the model and experimental 
values. 

The model curve for mean Stabilisation , with a 
step in temperature of -2°C, compares well with the observed 
stabilisation except below about 60 cm, where stabilisation 
is underestimated. The model curve indicates that the range 
of values of AT). in the six experiments was -0.5 to -9°C. 

Significant deceleration in the mean windspeed 
over the cool grass was observed in the average over 13 
wind-profile experiments. The deceleration was stronger 
than that predicted by the numerical model of Taylor (1970), 
but maximum deceleration occurred at approximately the 


height predicted by the model. 
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Sez Recommendations 

The following comments are intended to convey to 
possible future experimenters the lessons learned in the 
Ellerslie temperature-profile modification experiment. 

First and foremost, it is imperative that the 
System used to measure the mean temperature field be well 
tested and accurately calibrated prior to the experiment. 
Because the mean temperature differences are very small in 
comparison with turbulent fluctuations, it is impossible to 
Obtain definitive data with an improperly prepared network. 
Furthermore, even a satisfactory measurement system will not 
provide meaningful data unless the sampling procedure 
conforms with the requirements of stationarity and of 
sufficient number of data points. 

The Seeker dimension of the leading edge and 
modified surface should be made sufficiently large, or the 
downwind fetches sufficiently Small, so that at the furthest 
downstream mast there iS no possibility that the air 
arriving has not come across the leading edge. The site 
chosen must be large enough to allow this requirement to be 
met, and preferably large enough to allow all wind 
directions to be used with confidence that the upwind 
conditions are uniform, because the necessity of using only 
a subrange of wind directions may eliminate the possibility 
of experimentation On many otherwise satisfactory occasions. 

Instruments for profile measurements should be 


placed at the same heights on all masts involved, in order 
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to facilitate comparisons. For the same reason, the period 
of avelaging should be the same on all masts. 
Non-Simultaneous, spatially-separated profiles may be very 
difficult to compare, particularly if conditions are not 
very steady. 

In the Ellerslie experiment no Significant 
difference was observed between the average cooling on each 
of the downstream masts. In many cases , no. Significant 
difference was observed in windspeed between the two 
downstream masts. Because of this, it would seem logical to 
divide all available sensors between the upwind and a single 
downwind mast. This would be beneficial in optimising 
resolution (by which means it may be possible to determine 
the top of the IBL) and in simplifying the experimental 
layout. 

While ensuring that all necessary variables are 
measured, the experimenter should keep the experiment as 
Simple as possible, concentrating on performing the 
measurements with highest priority as accurately as 
possible. Every additional instrument used may demand 
pre-experimental work and planning and attention during the 
experiment which is better Spent on the primary 
measurements. 

If possible, upstream and downstream surface 
temperature should be measured, as the temperature step is 
usually required aS input for any model. 


Modification will be strongest under conditions of 
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strong insolation (and therefore strong evaporation). Such 
conditions are therefore preferable, as long as winds are 
not too light and variable. 

In conclusion, experiments of this nature are time 
consuming and may be expenSive. Consequently careful and 


unhurried planning and preparation is vital. 
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APPENDIX 


NOMENCLATURE 


Sensible heat exchange between surface and air, positive 


upward; (mW cm-?). 

Specific heat at constant pressure; (cal g-1 °K-1). 
Specific heat at constant volume; (cal g-! °K~-1). 
Evaporation rate, positive upward; (cm s~!1). 
Absolute humidity; (g cm-$3). 

Absolute humidity at saturation; (g cm-3). 


Surface absolute humidity upwind of discontinuity; (g 
cm—3). 


Surface absolute humidity downwind of discontinuity; (g 
Cis -)i- 


Change in surface absolute humidity (e,-e,)3; (g cm-3). 
Dimensionless function defined by equation (7.3.1). 


Soil heat flux density at the surface, positive upward; 
(mW cm). 


Acceleration of gravity; (cm s~2). 
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V7 
UE Incomplete gamma function ratio; (dimensionless). 
Ke Eddy diffusivity for water vapour; (cm? s~t). 
Ky Eddy diffusivity for heat; (cm2 s~!). 
Km Eddy diffusivity for momentum; (cm? s~!). 


K, Reference valve of diffusivity in power-law expression 
FOETALEfusivity so (eneesst)i. 


k von Karman's constant (0.4); also used for thermal 
conductivity. 


Pome 


Unit vector in vertical (z) direction. 


L Latent heat of evaporation of water, (585 cal g-!);also 
used for Monin-Obukhov stability length in Chapter 4. 


m Exponent in power law for wind profile; (dimensionless). 


n Exponent in power law for diffusivity profile; 
(dimensionless). 


Pp Pressure. 

q Specific humidity; (dimensionless). 
R Gas constant; (m2 s-2 °K~2). 

(RQ Rate of heat addition by radiation. 


R, Flux density of longwave radiation Cecieved at the 
surface, poSitive downwards; (mW cm-@). 


Flux density of shortwave radiation received at the 
surface, positive downwards; (mW cm-?). 
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Net radiation, positive downwards; (mW cm~2). 


Reflection coefficient of the surface for short wave 
radiation; (dimensionless). 


Temperature; (°K). 

Dewpoint temperature; (°K). 

Surface temperature upwind of discontinuity; (°K). 
Surface temperature downwind of discontinuity; (°K). 
Change in surface temperature, (T, -—T,) 3),(°K). 
Windspeed in x direction; (m s~!). 


Reference velocity in power law for velocity-profile; (m 
Satie 


FrLCcLOn wWELOCELY +. .(m (Ss; 2) 
Velocity vector; (m s~?). 
Windspeed in y direction; (m s~t). 
Windspeed in z direction; (m s~?). 


Horizontal coordinate perpendicular to leading edge; 
(cm). 


Horizontal coordinate parallel to leading edge; (cm). 
Vertical coordinate; (cm). 


Reference height; (cm). 
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Pp 


Constant of proportionality defined in Section 6.3. 


Constant of proportionality defined in Section 6.3. 


Height of IBL; (cm). 


Emissivity of surface for longwave radiation; 
(dimensionless). 


Equivalent temperature; (°K). 


Potential temperature; (°K). 


Rate of dissipation of kinetic energy into heat. 


Monin-Obukhov similarity function for momentun; 
(dimensionless). 


Monin-Obukhov similarity function for heat; 
(dimensionless). 


Transformation defined by equation (3.3.1). 


Incomplete gamma function; (dimensionless). 


Gamma function; (dimensionless). 


Density of air; (g cm-3). 


Density of liquid water; (g cm-3). 


Stefan-Boltzman constant; (cal s-! cm-2 °K—~*). 


Kinematic viscosity. 


(A prime denotes a turbulent fluctuation). 
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